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Abstract 

We study the inviscid multilayer Saint-Veuant (or shallow-water) system in the limit of small 
density contrast. We show that, under reasonable hyperbolicity conditions on the flow and a 
smallness assumption on the initial surface deformation, the system is well-posed on a large time 
interval, despite the singular limit. By studying the asymptotic limit, we provide a rigorous 
justification of the widely used rigid-lid and Boussinesq approximations for multilayered shallow 
water flows. The asymptotic behaviour is similar to that of the incompressible limit for Euler 
equations, in the sense that there exists a small initial layer in time for ill-prepared initial data, 
accounting for rapidly propagating “acoustic” waves (here, the so-called barotropic mode) which 
interact only weakly with the “incompressible” component (here, baroclinic). 
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1 Introduction 


1.1 Presentation of the models and the problem 

This work dedicated to the study of the so-called multilayer Saint-Venant system, which arises 
as an approximate model for the propagation of waves in the ocean or atmosphere, when density 
stratification cannot be neglected. We will refer to as free surface system the following first-order, 
quasilinear system of N(1 + d) coupled evolution equations: 


( 1 . 1 ) 


N 

dtCn + X! ^ = 0, 

< (n=l,...,7V) 

n 

5tu„ -k — - pj_i)VCi -b (u„-V)u„ = 0. 


Here, the unknowns C,n{t, x) and u„(f, x) represent respectively the deformation of the interface 
and the layer-averaged horizontal velocity in the n**' layer, at time t and horizontal position x € 
where d € {1,2}; see Figurej^ If d = 2, then we denote x = {x,y) and u„ = We denote 

by Pn > 0 the mass density of the homogeneous fluid in the layer, whereas g is the gravitational 
acceleration. Finally, 

/l„(t,x) Sn + Cn{t,X.) - C„+i(t,x) > 0 

is the depth of the layer. By convention, we set po = 0 (above the upper free surface is vacuum), 
and (N+iit,^) = 0 (the bottom is flat). We restrict ourselves to the setting of stable stratification, 
namely 

0 < Pi < • • • < Pat. 




Figure 1: Sketch of the domain and notations 


We may rescale the variables so as to replace the factor 


9_ 

Pn 


{Pr 


Pi-l) 


n 


In 


where 


def 'yi — 1 j def Pn 

rt = — - and 7 „ = —. 

1 - 7i Pn 


More precisely, we use the following nondimensionalization: 


X i x/A, C,n ^ Cn/ a, Sn Sn/ U, t i C^t/ A, llyj i Un/ Cq, 

where A is a characteristic horizontal length (say the wavelength of the flow), a is a characteristic 
vertical length (say the typical depth of one layer at rest), and cq a /(1 — 71)50 is a velocity. 

^Because we assume that the layers are all of comparable depth and the vertical stratification is balanced, in 
the sense that we fix m > 0 such that sup^^^^ ^, r“^} < m, then cq measures the typical velocity of 

propagation of the baroclinic modes; see Appendix Ib| 
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Although A does not appear in our system, it plays a very important role as the Saint-Venant 
system may be seen as the first-order asymptotic model obtained from the full multilayer water- 
wave system in the limit fi a} j}? —>■ 0; see [551 El] in the one-layer case and [T5| in the bilayer 
(albeit irrotational) setting. It may also be formally obtained using the hydrostatic and columnar 
motion assumptions; see [Ml El EH ESI . 

In this work, we ask 

Qn: What is the behaviour of the solutions to (|1.1[) in the limit 71 —)■ = 1? 


The first observation is that the velocity evolution equations become singular, as ri = 00 

since 70 = 0 by convention, so that even the existence of solutions on a non-trivial time interval is 
far from straightforward. 

At the linear level, it is known |3D] (see also Appendix that the flow may be decomposed 
into N modes, propagating as linear wave equations with distinct velocities. In our setting, the 
first mode, i.e. the barotropic mode, propagates much faster than the other, baroclinic modes, in 
the sense that the former is typically of size « ’^hile the latter are uniformly bounded when 

7 i —> I; hence the singularity. While such a decomposition is exact only in the linear setting, we 
show in this work that the flow behaves in a similar way in the weak density contrast limit even 
when strong nonlinearities are present, provided that the initial surface deformation is small, the 
depth of each layer remains positive and shear velocities are not too large. Roughly speaking, we 
show that one may then approximate the flow as the superposition of rapidly propagating acoustic 
waves and a non-singular “slow” mode with non-trivial dynamics. 

Let us be more precise. Asymptotically, the fast mode describes the propagation of the free 
surface, Ci, and total volume flux, namely IIw where w '=^ and n VA is the 

orthogonal projection onto irrotational vector fields. One has at first order, in dimension d = 2 and 
provided that Ci and IIw are initially balanced so that Ci « — 7 iw, the linear acoustic system 

(1.2) dtCi + V • (Hw) = 0; at(nw) + ^"=1 " VCi = 0. 

1-71 


The slow component contains all the baroclinic modes which interact strongly one with each other 
in the nonlinear setting. We show that it is asymptotically described by the rigid-lid model, which 
is obtained from the free-surface system by setting Ci = 0 in the mass conservation equations, and 
replacing riV^i with 'S/p in the velocity evolution equations (see |57||5|). In addition, we apply the 
so-called Boussinesq approximation (see e.g. [5D]) to the limit system, that is we set 7 „ = 1 in the 
velocity evolution equations while r„ remains fixed and non-trivial. This yields 


(1.3) 


N 

dtCn + X! ^ = 0 , 

i—n 

n 

dtUn S/p + (u„-V)u„ = 0, 

i-2 


(n = 1,..., A) 


with the same notation as before, except Ci(^) x) = 0 and hi{t, x) = (5i — f 2 {t, x). In particular, the 
first equation is not an evolution equation but a constraint (conservation of total mass), namely 


V • = 0. 

Physically speaking, p{t, x) is (up to a constant) the pressure at the flat, rigid lid. From a mathemat¬ 
ical viewpoint, p is the Lagrange multiplier associated with the above divergence-free constraint. It 
may be reconstructed from the knowledge of C 2 , ■ • •, Cwi Ui,..., u^r by solving the Poisson equation 



N 

Ap -f ^ V • 

71 — 1 



V)U„ -I- U„V • (hnUn) -I- hn 


i=2 


0 . 


We thus offer a rigorous justification (the formal justification is generally attributed to Armi |5|) 
of the widely used rigid-lid and Boussinesq approximations for free surface multilayer shallow-water 
flows with a small density contrast. 
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1.2 Main results 


Before stating our main results, let us recall for convenience the system of equations at stake: 


(1.4) 


N 


dtCn + X! ^ 


(n= 


1 1 

dtUn + — . VCi + —V r-iVCi + (u„-V)u„ = 0 . 

7n 1 - 7l 

where Si + Q — Ci+i with convention Ca^+i = 0 and ■ In the following, we fix 


parameters 5i,..., r 2 ,..., tat S ( 0 , oo); and denote 


def 


m = max {<5^,(5^ }. 

We can then reconstruct 7 i G (0,1) with = 1 — J2f=i+i ’’j (* = 2,..., N) and 71 = 1 — 
where g is the only parameter allowed to vary, and by assumption 

£» <C 1 . 


It is also convenient to denote (notice the g ^ prefactor; see Remark 1.3 below) 


U (g ^1X2, ■ ■ ■ ,Cn,u 7 , ... , 


SO that the control of U in Sobolev space (see Appendix [ a| for notations) yields, when d = 2, 

N N N 

Eh-lb + EKIb + EKI 


1 1 1 12 <^f 1 I > 12 

\^\h^ ~ 


n=2 


|2 

71—1 n—1 


Let us now state the main results of this work. 

Theorem 1.1 (Large time well-posedness). Let s > d/2 + 1 and U‘” G he such that 

(1.5) Vn G {1,..., iV}, inf /i“ > ho > 0, 

xGR*^ 

. dcf ’ * 

where h// = (5n + Cn” ~ Cn"+i) we recall the convention Cn+i = 0. 

One can set i2o,gQ^ = C{m,ho^, l'-^*”lArd that if li™ satisfies additionally 


( 1 . 6 ) 


VnG {2,..., IV}, sup |u“-u“_i| < Vq\ 

xGR'^ 


then there exists T > 0 and a unique U G {[0,T); strong solution to (1.4) and 

U = U“, and satisfying (1.5),(1.6) with /io/2 and r'o/2 for any t G [0,T). Moreover, one has 

T-i < xC(m,VMU“|^J and < |U“|^, x C(m, 

uniformly with respect to g & (0, gfi}- 

Theorem 1.2 (Strong convergence). Let d = 2, s > d/2 + 1, and U“ G as above. 

Then there exists T > 0 with T~^ < |U ™|^3 x C(m, /iq and, for g sufficiently small, 

• U G C°([0,T);iL®(IR‘^))'^(^+‘^) a unique strong solution to ( |1.4[ ) and U = U‘”; 


• (P^^) ) CV’ ..., u^'")^ G L°°{0,T; H‘^a unique strong solution to (1.3) 

with initial data 

Cn^L=o=C {n = 2,...,N) and = u“ - (n = 1,..., A^), 

where 5 total depth, w‘” J2n=ii^'>^ + Cn* ~ Cn'+i)'^™ convention 

= Cn +1 — 0 n VA“^V- the orthogonal projection onto irrotational vector fields. 
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^-l^ac^^ac g c°{[0,T);H^ 


c 


a unique strong solution to (1.2) with 
1 lt=o = C“ and 1,^0 = nw“. 


Moreover, one has for any 0 < s' < s, 

11^ ~ ® (£'—>' 0)i 

where U'^pp Cl^^> • ■ •. ..., 

Remark 1.3. As mentioned in the introduction, our hypotheses contain a smallness assumption 
on the initial deformation of the surface, namely Ci L^o = ^(p)- This assumptions is natural so as 
to balance the contributions in the (preserved in time) energy: 

f + '^rn\Cn\^ + '^-fnhn\Un\^. 

^ n=2 n=l 


Without this assumption, the flow possesses a strongly nonlinear barotropic component, and energy 
methods yield a well-posedness theory over a small time-domain, t G [0,T), T~^ = 0{g~^); see 
Proposition \2. 1\ and Remark 2.2. below. On this timescale, the baroclinic component do not evolve, 
so that all the dynamics is described by the barotropic component (asymptotically as g ^ 0). 

Remark 1.4. The requirement g G (0, go) does not lose in generality in Theorem [HI- the case of 
non-small g follows from the standard well-posedness theory of quasilinear systems, as proved by 
Monjarret in m and stated in Proposition\2. 1\ below. 


Remark 1.5. Our proof does not rely on, but rather provides, the existence and uniqueness of strong 


solutions of the limit (rigid-lid) system. In that respect, one may see the free-surface system (1.4) 
as a penalized model for ( |1.3| ) relaxing the rigid-lid constraint. Sharper well-posedness results for 
the rigid-lid system in the two-layer case and without the Boussinesq approximation are provided 
in \^. 


Remark 1.6. Theorem 1.2 is restricted to d = 2 because we use dispersive decay estimates on 


rapidly propagating acoustic waves in order to control nonlinear coupling effects between the fast 
and slow modes. In the case of dimension d = 1, and provided that the initial data is sufficiently 
localized in space, we justified in m a similar mode decomposition of the flow, by making use of 
the different spatial support of each mode after small time. Proposition 4-4 therein, together with 
Proposition \3.^ in the present work, offer a convergence between the exact and the approximate 
solution with rate 0{g). The same convergence rate holds in the case of dimension d = 2 and 
well-prepared initial data, in the sense of Proposition \4^ 

In the following, for the sake of simplicity, we limit our study to the case of dimension d = 2, 
although we find it more telling to keep the notation d. The proof of Theorem \n\ is easily adapted 
to the case of dimension d = \. 


Remark 1.7. One could add, without any additional difficulty, a uniformly bounded and order-zero 
term to the system, so as to take into account for instance the Coriolis force, atmospheric pressure 
variations, or bottom topography. Notice however that these terms should be of size 0{g) in the 
evolution equation for Ci; in particular, only small topography may be dealt with using directly our 
strategy. Similarly, except in the one-layer case where the the component due to Coriolis effect is 
an antisymmetric perturbation of a symmetric system, one cannot allow a rapid rotation such as 
in the quasi-geostrophic regime, which would correspond here to Ro « g where Ro is the Rossby 
number; see II3 [Si- 

Remark 1.8. Our results are valid for arbitrary N, but not uniformly. In particular, we cannot 
control the dependence of go, vq as N grows. Thus our strategy cannot be adapted to study the system 
in the limit N ^ oo, corresponding to the physically relevant situation of continuous stratification. 
A similar shortcoming was already noticed and discussed by Ripa in \37(. 
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1.3 Discussion and strategy 

A well-posedness result on system ( |1.4[ ) is stated and proved in m- It follows from a standard 
analysis on quasilinear systems since a symbolic symmetrizer may be exhibited (see Sectionj^below). 
However, due to the presence of singular components in the system, the a priori maximal time of 
existence of the solutions may be bounded from below only as T^ax ^ Q using this method. Such a 
result is unsuitable for our purpose as the time interval shrinks to zero in the considered limit, and 
is inconsistent with oceanographic observations of large amplitude internal waves propagating over 
long distances; see e.g. [33] and references therein. 

In order to go beyond this analysis and provide a time of existence of solutions uniformly 
bounded from below with respect to g small, we need to take advantage of some additional structural 
properties of the system. This structure is put to light by a suitable change of variable, which we 
describe below. 

Let us introduce the shear velocities, v^, and the total horizontal momentum, w, as follows: 


N 


(1.7) 


w 


^ Kun and Vi yiU* - 7 ^_iUi_i (z = 2,..., iV); 


SO that, conversely, for any n = 1,..., A/", 

N 


( 1 . 8 ) 


u„ = 


7n^ 


” “ -U 




with 


ifj<n. 




^hi otherwise. 




One may rewrite system (1.4) using the new variables as follows: 
dt{g~^C,i) + -V • w = 0, 


N 


(1.9) 


dtC,n + X! V • (hiVi,) = 0, 

i—n 

dtVn + TnVCn +7n(Un ’ V)u„ - 7 „_i(u„_i • V)u„_i = 0, 


(n = 2,...,7V) 


N 


9tw+ I ^7^- 


7i 


N / N 


N 


1 - 7i 




i=2 \j-i 


i=l 


where V • (hiUi 0 Ui) ^ ‘ 


and Ui = (uf,uf)^ are meant as the expressions 


' d^(hilufl^)+dy(hiufuf) 

. dx (^iuf uf )+dy (hi juf ) 

in terms of (Ci,. •. Cn,'^ 2 , ■ ■ ■, vv, w)^ given in ( |1.8| ). 

The above change of variables may be seen as an approximate normal form allowing to decouple 
the slow and fast components of the flow. Indeed, since 


N 


N 




i=i 


one sees immediately that the singular terms appear only as linear components on the evolution 
equations for s-nd w —or more precisely Hw— and involve only p“^Ci Hw. In other 

words, the leading-order terms form a system of rapidly propagating acoustic waves in 

dtig~^Ci) -b -V • w = 0, 

Q 

5zw + i(Ef=i<5,)v(p-^Ci) = o(i). 

The remainder contains quasilinear components depending on both the fast (p^^CijIIw) and slow 
(C 2 , ■ • •, Cn) V 2 ,..., v„, (Id — n)w) variables, so we need to consider the full system of equations (1.9) 
in order to obtain the desired uniform energy estimates. 
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Let us now, for the sake of simplicity, restrict our discussion to the case of dimension d = \. 


From the above, we may rewrite the system (1.9) as 

( 1 . 10 ) 


atV + = dtM + -L9,rV + CMa^rV = 0 

Q Q 


where V {g ■ ■ Xn ■ ■ ■ ,vn , so that 9tV + ^Li9a;V = 0 represents the above 

acoustic wave system for the fast variables, while ker(L) = span(C 2 ,..., Cat, V 2 , ■ ■ ■, vn); and C[V] 
contains lower-order (in terms of g) and coupling terms. 

We shall make use of the fact that one can construct a “good” symmetrizer of the system under 


the form (1.10), namely we exhibit real, positive-definite matrices T[V] such that T[V] = (T[V])^ 
and T[V]B[V] = (T[V]B[V])^, and satisfying the decomposition 

T[V] = T(o) + T(i)[(ld - nf)V](ld - Hf) -f Oig), 

where (Id — Hf) is the orthogonal projection onto ker(L), the slow variables. Indeed, one obtains an 
energy estimate by taking the inner-product of the equation with T[V]V, which only requires to 
estimate 

|5,(iT[V]B[V])|^„„ + |a*(T[V])|^„„. 

Using that L and T(g) are constant operators and that (Id — FlfjcltV = —(Id — nf)C[V](92,V = 0{1), 
we see that the above are estimated uniformly with respect to g small; thus we have a uniform 
control of the norm. The corresponding estimate with s > d/2 + 1 does not bring additional 

difficulties, using that L commutes with the Fourier multiplier (1 — Aj'^/^ld. 


There remains to understand why such symmetrizer exists for our system (1.10). One could 
check, after tedious calculations, that the explicit one provided by Monjarret in m (after apply¬ 
ing the congruent transformation associated with the change of variables) satisfies the necessary 
hypotheses, but we offer in Appendi x [B] an alternative and more robust construction. We show 
that, provided that V satisfies (1.5),(1.6), then ^B[V] has 2N real and distinct eigenvalues. Two 

of them are asymptotically equivalent to A± « ± ^ Sj as > 0 while the other ones are 

uniformly bounded with respect to g small. The spectral projection corresponding to the former 
converge towards the projections onto the eigenspaces corresponding the two non-trivial eigenvalues 
of L. Using the scale separation between the eigenvalues, one shows that the spectral projections 
corresponding to the latter are uniformly bounded, and that they converge as g —> 0 towards 
independent, rank-one projections onto subspaces of ker(L). Our symmetrizer is then, classically, 

2N 


T[V] ^(P, [V])^P,[V] 


1=1 

where Pj[V] is the spectral projection onto the eigenspace of B[V]. That T[V] enjoys the de¬ 
sired properties follows, using standard perturbation theory |22j . from the fact that L is constant, 
ker(L)0_Lran(L) = and the strong scale separation between |A±| > 1/g and uniformly bounded 
eigenvalues. 

An additional difficulty arises in the situation of horizontal dimension d = 2, due to the fact that 
the symmetrizer of the system —which is constructed from the symmetrizer in dimension d = 1 and 
a rotational invariance property— is only a symbolic symmetrizer, as opposed to symmetrizers in the 
sense of Friedrichs. Thus we rely on para-differential calculus, but extra care must be given to “lower 
order terms” in the sense of regularity, which may effectively hurt our energy estimates if they are not 
uniformly bounded with respect to g. As a matter of fact, we use that one can construct an explicit 
operator (defined as a Fourier multiplier) which symmetrizes the linear, singular contributions of 
the system, and use para-differential calculus only on the next order components in terms of g. 

Given the uniform (with respect to g small) energy estimates, the large time well-posedness 
(Theorem |l.l[ ) follows from the standard theory for quasilinear hyperbolic systems. The convergence 
results (Theorem 1.2 as well as additional assertions in Section]^ proceed from rather standard 
techniques in the study of singular systems; see references below. 
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1.4 Related earlier results 

In |16j . the author studied the so-called inviscid bilayer Saint-Venant (or shallow water) system in 
the limit of small density contrast. The change of variables allowing for uniform energy estimates 
was exhibited therein, and convergence towards a solution of the rigid-lid limit, as well as a second- 
order approximation, was deduced in the case of well-prepared initial data. This work is therefore 
an extension of these results to the situation of (horizontal) dimension d = 2, ill-prepared initial 
data as well as arbitrary number of layers. 

As already noticed in the aforementioned work, our problem has many similarities with the 
(two-dimensional) incompressible limit for Euler equations, as studied initially in As 

a matter of fact, if we consider only one layer of fluid, then our problem corresponds exactly to 
a special case of the isentropic incompressible limit, and we recover the results of Ukai [22] and 
Asano [3|. We will not detail the very rich history of results concerning this problem (we let the 
reader to [1911281 |T| for comprehensive reviews) but rather aim at pointing out similarities and 
differences of our situation. 

Let us first recall the two-dimensional isentropic Euler equations for inviscid, barotropic fluids: 


f dtp + V ■ (pv) = 0, 

\ 9t(pv) -I- ^VP -I- V • (pv (g) v) = 0, 


where P = P{p) is a given pressure law, p > 0 is the density, v the velocity, and e the dimensionless 
Mach number. As claimed above, one recognizes exactly (1.91 in the one layer setting {N = 1), by 
setting P{p) oc p^, and identifying 


^1 — ^1 + Cl 


pv 


and 


p. 


Of course, the difficulty in our case is that, as one considers additional layers of fluids, these equations 
are coupled with additional equations on additional unknowns, so as to produce a full quasilinear 
system. Since these additional equations are non-singular with respect to the small parameter, it 
is tempting to compare our situation with the incompressible limit for the non-isentropic Euler 
equations, where (1.11) is coupled with an additional evolution equation for the entropy, S\ 

dtS+{vW)S = 0 


and P = P{p, S). 

Our situation, however, is quite different. This can be seen from the fact that, contrarily to 
the non-isentropic Euler equations, the linearized system is balanced, in the sense that a small 
perturbation of the “slow” component of the reference state induces only a small deviation for the 
solution. In other words, using the notation of the above discussion, the symmetrizer of the non- 
isentropic Euler equations does not satisfy T[V] = T(q) -|-T(i)[(ld — nt)V](ld — Hf) +0{g) but only 
T[V] = T ( 0 ) +T(i) [(Id— nf)V] -|-C>(p); see discussion in |3D|. This additional property in our situation 
allows in particular to straightforwardly deduce energy estimates from the corresponding 
energy estimate; and to obtain the strong convergence result of Theorem 1 1. 2| simply from dispersive 
estimates on the “acoustic” component of the flow, as originally carried out by Ukai [42] and Asano |3] 
in the isentropic case. 

In order to deal with this situation, Metivier and Schochet (SD) (see also [7]) rely on the fact that 
their system enjoys a diagonal block structure and that the symmetrizer commutes exactly with 
the singular operator, denoted ^ \-dx in the above discussion. Roughly speaking, this means that we 
may control the compressible and isentropic component of the flow independently of the acoustic 
component by simply projecting onto ker(L). Since such assumptions are only approximately satis¬ 
fied in our situation, our system is rather related to the “p-balanced” (and not “p-diagonal”) systems 
studied by Klainerman and Majda |23j . although we do not restrict ourselves to well-prepared initial 
data. 

In this spirit, our proofs rely as little as possible on explicit calculations, thus we expect that 
the general strategy may be successfully applied to other situations and other frameworks, such as 
the ones presented in |3S| . This is why we use mostly in the following the terminology of “fast vs 
slow” mode/component instead of “barotropic vs baroclinic” which is more relevant to our initial 
oceanographic motivation (see [ID]); or “acoustic vs incompressible” associated with Euler equations. 
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1.5 Outline of the paper 

In Section]^ we recall that our quasilinear system admits an explicit (symbolic) symmetrizer, which 
yields immediately a well-posedness theory for (1.4), for any fixed p > 0. 

In Section |3.1l we exhibit the structural properties enjoyed by our system, after the change 


of variables ( 1.7),( TT^, and which allow the uniform (with respect to g small) energy estimates 
provided in Section |0| (Proposition [3^ and Section 3.3 (Proposition |3.8[ ). 

We deduce in Section the large-time well-posedness result of Theorem o Additionally, we 
show in Proposition |4.2| that an assumption of well-prepared initial data is propagated by the flow 
for positive times. 

Sectionj^is dedicated to convergence results. We first state in Proposition |5.1| a weak convergence 
result for the solutions of the free-surface system as p —>■ 0. As in the incompressible limit for Euler 
equations, the convergence cannot be strong uniformly in time, due to the rapidly propagating fast 
mode. However, we show in Proposition |5 . 2| that this small initial layer in time vanishes in the case 
of well-prepared initial data, and then characterize the defect for general initial data in Section [5.3[ 
yielding Theorem ] 1.2 [ 

Appendix contains a description of some notations used throughout the text, as well as a 
short review of standard results concerning product and commutator estimates in Sobolev spaces 
(Section [A.2[ ) and Bony’s paradifferential calculus (Section [A.3[ ). 

Finally, Appendix is dedicated to some results on the eigenstructure of our system, which are 
used in Section IXTl 


2 Standard well-posedness theory 


In this section, we fix p > 0 and construct an explicit (symbolic) symmetrizer of (1.4). This offers a 


well-posedness result similar to Theorem ] 1.1 [ although non-uniformly with respect to g small. This 
analysis has been provided by Monjarret in [31] : we recall it here for the sake of completeness and 
because the objects defined therein will be of later use. 

Let us first rewrite the free-surface system (]1.4|) with a matricial, compact formulation. Provided 


that U ‘'A* (p-^Ci, C 2 ,..., c^, , ulY e CO([0, T); H\W ^))^with s > d/2 + 1 


one can rewrite (1.4) equivalently as 


( 2 . 1 ) 

with 


dt\} + -A^[U]9^U -b -A2'[U]9,,U = 0, 
g g 


( 2 . 2 ) 


-A’’’[U] 

0 



H(C) I 

On \ 

:a: 

[□(«“)] 

Oat 

V On 

! 0^ ! 

D(«7 / 


^A^U] 

0 



Here and thereafter, we heavily make use of the block structure of (l-l-(i)-by-Af (l-bd) matrices. We 
denote by Oat the N-hy-N matrix with only zero entries, and for u G D(u) = diag(ui,..., mat). 
Moreover, M, H are upper-triangular and R is lower-triangular and are defined by 


M(u)„,i 


^(Ui - Uj_i) 
Ui U-i—i 
Ui 

0 


if 1 = n < z, 
if 2 < n < z, 
if i = n, 
a i < n, 



1 71 
Q In 

Vj 

7rr 

0 


if n = 1, 
if i > n > 2, 
if i < n, 

if j = 1, 
if n > i > 2, 
if n < i. 
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Proposition 2.1 f|31j. Theorem 2.8). Let s > d/2 + 1 and U“ £ be such that 

(2.3) Vn e {1,.. . ,7V}, inf = inf {Sn + (n - Cn+i) > ho > 0, 

xGR"^ xGR'^ 

where we recall that Cn+i = 0 by convention. One can set vq = ^ such that if U™ 

satisfies additionally 

(2.4) yn£{2,...,N}, sup |< - <_i| + |<-<_i| < 


then there exists a unique Tmax > 0 and U € C°([0, Tmax) ; maximal solution to ( | 2 . 1 [ ) 

and U = U“. 

Moreover, if < oo, then |U|j^i ^(f) —>■ oo fit f T^ax^ or one of the hyperbolicity condi¬ 


tions (2.31,(2.4) ceases to be true. 


Remark 2.2. Naively following the above strategy and keeping track of the dependency of constants 
with respect to the parameter g would yield a disappointing lower bound on the maximal time of 
existence, namely Tff/y. < , even when the initial surface deformation is assumed small as in 

Theorem o It is the main result of this work that, in this case, the well-posedness theory and 
uniform energy estimates can be extended to a non-shrinking time domain as g ^ 0. 

The conditions (2.3),(2.4) ensure that the symmetrizer, defined in drzt below, is coercive. It 
is therefore a sufficient condition for hyperbolicity. Except in very specific cases, one has very few 
information on the domain of hyperbolicity; see discussion in Appendix\^ 

The fact that (2.4) requires a control on the shear velocities, u„ — u„_i rather than on the veloc¬ 


ities themselves is allowed by some freedom in the choice of the symmetrizer. Notice in particular 
that the Hessian of the conserved energy yields a natural symmetrizer for our system of conservation 
laws in the case of irrotational flows, but it does not enjoy the desired property. 

The ability to construct symmetrizers depending strongly on the shear velocities but only weakly 
on a background velocity (or on the total volume flux) is also essential for us to obtain results outside 
of the scope of well-prepared initial data; see m for an analysis where this property is not used. 


Proof. Let us introduce the symbol of (2.1): 


A[U,e] ‘^= rA“[U] +eA*'[U]. 

An important property of our system is that is satisfies rotational invariance [?n Section 1.2]. More 
precisely, one can easily check that 


(2.5) 


ve = (r,e)^ e R" \ {oj, A[u,e] = Q(o-'A"[Q(e)u]Q(c)ici, 


where 


( 2 . 6 ) 


Q (0 


( ICI'w 

! 0^ 1 

On \ 

Ojv 

1 


V Ojv 

!-eiiv! 

C\n / 


where I AT is the 7V-by-A^ identity matrix and 1^1 = + Obviously, Q(^) is homogeneous 

of degree 0 in with entries in \ {0}) and is orthogonal: Q(^)“^ = Q(C)^- 

This allows to construct a (symbolic) symmetrizer of the system from a (Friedrichs) symmetrizer 
of alone. More precisely, define 


(2.7) 


S^[U] 


^__D(r)__l __V__! \ 

_ _ _L _ _ rojy) P _ _0^_ _' 

y Oat j Oat j D(7)D(/i) y 
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with D(f) = diag( 7 i,r 2 ,... ,rAr) and L = -D[ 7 (u“+ itr^)]AD[eg] where D(eg) = diag(gi, 1,..., 1), 


A 

/-I 1 

(0)\ 

, A-i = 

(-1 

(-1)\ 


1(0) 

1 

-l) 


\ (0) 

-1 y 


and the parameter may be chosen freely in span(rt“). 

That S^[U]A“[U] is symmetric is easily checked once we notice the following identities in (2.1 1 : 
M(u) = D(eg)-iA^iD(u)AD(eg), H = -D{eg)-^A^^D{h) and R = -D( 7 -i)(A-i)^D(f)D(eJf^. 

It follows that S[U,^] Q(^)“^S^[Q(5)U]Q(^) is a symbolic symmetrizer of (2.1), since 

S[U,^]A[U,e] = Q(e)^S"[Q(OU]A"[Q(OU]Q(Oie| 


is obviously symmetric. Finally, one may choose = —uf for instance, so that as soon as U 
satisfies ( 2 .3),( 2 .4) with vq sufficiently small, then S^[Q(^)U] is strictly diagonally dominant with 
positive diagonal entries, and therefore definite positive, uniformly for any ^ G \ {0}. 

Proposition |2.1| is now a direct consequence of the standard theory of first-order hyperbolic 
quasilinear systems; see [29] for instance. □ 


3 Uniform energy estimates 

In this section, we establish uniform (with respect to g small) energy estimates, which are the 
essential ingredients in the proof of our main results. We first exhibit in Section [3T| some properties 
of the system obtained after the change of variables (fFTt-drst, and which allow the energy 


estimate of Proposition 3.5 (Section 3.2) and in turn the 77® energy estimate of Proposition 3.8 
(Section |3.3| ). 

3.1 A new formulation 

In what follows, we fix Aq > 0 and always assume that 

(3.1) Vu G {1, . . . , N}, hn = Sn + Cn — Cn-|-1 > Aq > 0 

(recall the convention: (n+i = 0). More precisely, we work with Yho C defined by 

V/io = |(£'“^Ci A2,..., Civ,*,. ■. ,*)^ G ]S^N{i+d) pjj ) holdsj. 


Consequently, the change of variables (1.7)-(1.8) define self-homeomorphisms between 


u — (e ^Ci) C 2 , • ■ •, Cn ,^ 


, . . . , tijY? 5 • • • ? ^j\[ 


GY^ 


ho 


and 


V (g \i,C 2 , ■ ■ ■ ,Cn,V 2 , ■. ■ ,v%,w^,vl^,... G Vfc 


and we denote 


'^ho '^ho 

V U F(V) 


and F ^ : 


V/io Yho 

U V=^F-i(U) 


Consider the Jacobian matrix associated to 7^ ^ 


Itv I Otv I Oat \ 

(:(u^)| 0 ^ |a,(c )7 


(3.2) 
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where 


/- 


A/.(C) = 


7i 72 


_( 0 )_ 

V hi 


( 0 )\ 


-7Af-i In 
Hm / 


and C(it) = 


( 0 ) 


y gui U2 — Ui ■ ■ ■ UN — UN-1 J 


From the inverse function theorem, one has 

(3.3) J^[V] = (J^-yF(V)])-i 

and 


\n 

~\-lr 


(J^- [U])-1 = I 

-A^icr^ciuy) 

where (as is easier seen directly from (| 1 . 8 |)), 




On 

On 


On 

On 

~A~icA 


E N — 1 l 

^=ll^ h 


with i X)Li 7 , ^h. 


if j = AT + 1 , 
if 2 < j < n, 


2=1 H 

“Z)i=j 7 i~^^i otherwise. 


Applying the change of variables U = F'(V) in (2.1) yields, for sufficiently regular functions (see 


Lemma 4.1 below), 
(3.4) 


which we can identify with ( |1.9| ) , and where [V], B*' [V] are explicitly given in terms of [V], [V] 

(displayed in \2.2\) and J^[V] through 


dtV + - B" [V]a, V + - B^ [V]9, V = 0, 

Q Q 


(3.5) B“[V] = (J^[V])-iA^[F(V)]J^[V] and B^^y] = (J^[V])-iA2'[i^(V)]J^[V]. 

Finally, we denote by Hf, Hf the orthogonal projection onto the “fast” variables, namely 


(3.6) 


def 


□f = 


D(ei) j Oat j Oat \ 
Ojv ~["D(eAr)~] Oat 
On I Oat jD(eAr)/ 



/ 

D(e0 

Oat 

On \ 



Oat 

|D(eAr)”[ 

Ojv 


1 

Oat 

I 0^ j 

Ojv / 


with D(ei) diag(l, 0,..., 0) and D(e 7 v) diag(0,..., 0,1). 
The following result is now straightforward. 


Lemma 3.1. Let U = A(V) € V/jq, and recall the definition o/Q(^) in (2.6). Then 
(3.7) V^gK‘'\{ 0 }, Q(C)F(V) = F(Q(e)V) and Q(OJ^[V] = J^[Q(OV]Q(e). 

Moreover, J^[V], (J^[V])“^ : Yho —>■ A1 Af(i+(i)(I^) are well-defined and smooth and satisfy 


(3.8a) 

(3.8b) 

(3.8c) 


|7"(V)| <C(m,VMv|)|V|, |T-i(U)| <C(m,VMu|)|U|; 


|J^[V]|| <C(m,hoMv|) 


|(J^[V])-i|| <C(m,/ioMV|); 


|J^[V] - J^[(ld - nf)V]|| < g C(m, |V|)|V|. 


Proof. Only the last estimate requires an explanation. Remark that the first variable of V con¬ 
tributes to Ah{C)~^ only through hi = — ^2 + g^- Thanks to the g prefactor, one deduces 

that 

||Aa(C)-' - AA((ld - nf)C)-i < p C{m,h^\ |C|)|Ci|, 
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where we denoted (Id — Hf)^ (0, C2, ■ • ■ Xn)^■ Similarly, by (1.7), one has 

Ui - Ui-i =Vi + {l- + (1 - 

and, by definition, 1 — 7i = ^?^ 'l2f=i+i so that 

||C(u(V)) - C(u((ld - nf)V))|| < e C(m, |v|)|v|. 


where u(V) represents the velocity variables of ^(V), as given by (1.8). The result is now clear. □ 

Lemma 3.2. The functions B“, : Yha —>■ AfAr(i+d)(®) o.f's well-defined and smooth and satisfy 

the rotational invariance 

(3.9) VC G M-" \ {0}, B[V,C] =' rB"[V] + eB^[V] = Q(C)-iB"[Q(C)V]Q(C)|C|; 

as well as the following estimates: 

(3.10a) ||B^[V]|| < C(m,/io\|V|), ||(ld- 

(3.10b) ||B"[Vi] - B"[V2]|| < Q C{hf,- 


B"^[V] 

< Q C(m,/io 

,|Vi|, 

|V 2 |)|Vi-V 2 | 


Proof. Identity ( |3.9| ) is deduced from ( |3.5| ), (2.5) and ( |3.7[ ): 

C“B“[V] +C^B?'[V] = (J^[V])-i(C^A“[F(V)] +C^A2'[J^(V)])J^[V] 

= (J^[V])-iQ(C)-^A-[Q(C)F(V)]Q(C)J^[V]|C| 

= Q(C)-1(J^[Q(C)v])-ia-[f(Q(c)v)]j^[Q(C)v]Q(c)|CI 
= Q(C)-'B“[Q(C)V]Q(C)|C|. 


Estimates ( |3.10| ) may be deduced from identity (3.5) and the explicit expressions of A®[F(V)] and 
J^[V], (J^[V])“^ given in (2.2),(3.3),(3.2). They are also apparent when identifying (3.4) with (1.9). 
In particular, one immediately sees that the evolution equations for C2, • • • j ^2,..., vjv are non¬ 
singular (for g small), and the singular term on the evolution equation for w"^ involves only dyfi, 
so that 

||i(ld-nnB"[V]|| <C(m,/ioMv|) 


and 


|-(ld - n?^)B"[Vi] - -(Id - nnB"[V2]|| < C(V\m, |Vi|, |V2|)|Vi - V2I. 
Q Q 


The only singular terms (for q small) arise from the first and last equations, which read 

-V(£<“^Ci) = rl\\/]dXJ + r5^[V]5yV, 


N 


dt{e ^Ci) + Q • w = 0 and dtw + ( 7^ ^hj 

U'=i 


1 - 7i 


where rf,rf are smooth and enjoy the same estimates as -(Id — nf)B“ above. We now only need 

^ .2 


to remark that = liQ ^ (by definition) and, since 1 — 7^ = 


-7i 

N 


N N ^ N 

l-7j 


H H H 7i) + 072(V), 

i=i i=i i=i i=i 


where, again, r 2 is smooth and uniformly estimated as above. Lemma[3.2|is now straightforward. □ 


Finally, let us provide some estimates on the symbolic symmetrizer of the system. 
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Lemma 3.3. There exists : Yhg A^jv(i+d)(I^) ® smooth function such that 

VVgV?,„, (T"^[V])T =T"^[V] and (T“[V]B^[V])T = T^[V]B^[V]. 

Moreover, one can set gQ^,n = C{m,hQ^, |V|) > 0 such that if g G (0, po) and \/ ,\/i ,\/2 satisfy 


(3.11) 

VnG{2,...,iV}, 

sup \vf\ + 
xeR'^ 

then one 

has the following estimates: 


(3.12a) 

cold<T"[V], ||T"[V]|| <6-0, 

with Cq^ = 

(3.12b) 

||n[Vi]-T"[V2]|| + -||(n[Vi] 
Q 

-n[V2])nf"|| 

(3.12c) 

||T“[Vi]-T"[V 2]|| <C(V\m, 

|ViMV2|)(|(l 


Proof. One could define the symmetrizer as T®[V] = (J^[V])^S^[F(V)]J'^[V], where S^[i^(V)] is the 
symmetrizer associated with A'^[F(V)] and has been displayed in (2.71, and check the properties 
directly on this explicit symmetrizer. Some of the estimates, however, rely on delicate cancellations 
which have no obvious explanation using this method. This is why we find it more instructive to 
construct our symmetrizer using the spectral properties of our system, whose study we postpone 
to Appendixfor the sake of readability. It is proved in Lemmas B.l and B.3 that B“[V] has only 
real and semisimple eigenvalues: 

1 ^ 

-B"[V] = ^ Mn[V]P„[V] + /r_„[V]P_„[V] + <[V]P° [V], 


where /Lt±„ [V], [V] G K, and P±n[V], P°[V] are rank-one spectral projections. We now define 


N 


(3.13) 


T-[V] = ^(P„[V])Tp„[V] + (P_„[V])Tp_„[V] + (P°[V])TpO [V]. 


n—1 


Indeed, T®[V] is obviously symmetric, and so is 


1 ^ 

-T"[V]B"[V] = ^A^„[V](P„[V])Tp„[V] -PM_„[V](P_„[V])Tp_„[V] +<[V](P°[V])^P0[V], 

n 


n—1 


and one has for any W G 


N 


(T"[V]W,W) = ^|P„[V]W| + |P_„[V]W| + |P°[V]W| > 


1 


{N{i + d)y 


IwP. 


The upper bound in (3.12a) follows from 


N 


^||P„[V]|| + ||P-„[V]|| + ||pO[V]|| < C(m, |V|), 


n—1 


which is given by (3.8bI in Lemma 3.1 with (B.l) in Lemma B.l and (B.2c) in Lemma B.2 
with (B.3c)-(B.3e) in Lemma B.3 This proves (3.12a|. 

Similarly, (3.12b) and (3.12c) are easily deduced from the explicit expression for P° in (B.l), 
and the estimates (B.3c),(B.3d),(B.3e) as well as (B.2e). □ 


We conclude this section by collecting the above information on the symbols of our operators: 
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Corollary 3.4. Let V € Who satisfying (|3.11[) and ^ \ {0}, define 


B[V,C] = rB"[V]+eB*^[V] = Q(e)-iB"[Q(0V]Q(0iei; 

T[V,e] = Q( 0 -'T"[Q( 0 V]Q(C) and Pf(C) = 

Then one has (T[V,^])^ = T[V,f], (T[V,^]B[V,^])^ = T[V,^]B[V, and the following estimates: 


(3.14a) 

(3.14b) 


|B[V,e]|| <Co Id and 
cold<T[V,d 


||(id-Pf(d)B[v,d|| <eCo Id; 
and ||T[V,d|| < Co; 


with Cg ^ 


Cn * = iV(l + d), Co = C{m, /ig , |V|); and for any Vi, V2 G Yfio satisfying ( |3.11| ), one has 

(3.14c) ||B|Vi,d-B[V2,d|| <eCi IV1 -V 2I Id; 

(3.14d) ||T|Vi,dB|Vi,d-T|V2,dB[V2,d <eCi Vi - V 2 |d; 


(3.14e) 

(3.14f) 


||T|Vi,d - T[V 2 ,d|| < Cl (|(ld - nf)(Vi - V 2 )| + e|Vi - V 2 I) ; 
|(T|Vi,d-T|V2,d)Pf(0ii <eCi IV 1 -V 2 I . 


mt/iCi =C(V\m, |Vi|,|V2|). 

3.2 energy estimate 

The following proposition shows that, thanks to the structure of the system exhibited in the previous 
section, one is able to control the (L^) energy of solutions, uniformly on a time interval independent 
of g small. This is the key ingredient in the proof of our main results. 

Proposition 3.5. Let V G 1T/^((0,T) X satisfying ( |3.1| ),( [04] ) with ho^v > 0, and 

W G C°([0,r);L2(K‘^))^(i+‘^), R G ^^(0, T; L2(Kd))Ar(i+d) ^^g^ that for any t G [0,r), one has 

(3.15) 

One can set ,vo = C(m, /ig ^ 


dt\N + iB"^[V]a,,W + -^&[Y]dyy\l = R. 


'-0 > ir llL<»((o.T)xB‘i)) if Q ^ (0, £>0) and v > vo, then 


(3.16) 


WU(0<Co(0)e' 


Cl y t 


W 


t=0 1^2 


„Ci V (t-t') 


with Co{t) = C(m, /ig |V| 


wi 


Ivlll llvl 




(t)), Cl = C(m, V' 

ell^tV 


V 


L“((0,T)xI 


L“((0,T)xR‘*)' 

+ ||(id-nf)5tV| 


Co(t')|R|i2(i') dt', 


0 and 


L“((0,T)xI 


This section is dedicated to the proof of this result. The main ingredients are the properties of the 
symbol of system (3.15) as well as its symmetrizer, collected in Corollary 3.4 Energy estimates for 
such symmetrizable systems can be obtained thanks to Bony’s paradifferential calculus associated 
with these symbols; see [29]. We shall however be cautious as paradifferential calculus typically 
provides estimates “up to lower order operators”: while this is sufficient for regularity aspects, 
this could induce order-zero but large (namely non uniformly bounded with respect to g small) 
remainder terms, preventing the desired uniform energy control stated in (3.16). 

This is why we decompose the symbols into a first order contribution which admits a natural 
quantization as a Fourier multiplier, whereas only second-order contributions will be paradifferen- 
tialized. Let us be more specific. Define 

6B{t,K,0 "=^'x(0<5B"|Q(0V(t,x)]|C| "='x(C)(B"|Q(0V(t,a:)] - B"|0])|C| 

where x is a smooth non-negative cut-off function {xiO = 0 for |^| < 1/2 and xiO = 1 for ICI ^ 1); 
and let TiSB be the associated paradifferential operator (see Definition A.2). Similarly, define 


<5T(t,x,C) =' x(a<5T“[Q(0V(t,x)] x(0(T"|Q(e)V(t,x)] - T"|0]) 


def 
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and Tst the associated paradifferential operator; and 


(3.17) 


s{t) (Q(n))- 


T"[0]+2('r5T + r/T) + AA-i)Q(i7), 


where ^ = (Id + \D\^) ^/^((Id — Hf) + pFIj) and A > 0 will be determined later on. 


We claim in Lemma [3)6| below, that the properties on T[V, given in Corollary 3.4 are sufficient 
to show that 5(t) is a uniformly bounded and coercive operator, and show a precise estimate on 
its time derivative, S'(t). In Lemma 3.7 we then rewrite (3.15) as a symmetric, paradifferential 
equation, from which the energy estimate (3.16) is easily deduced. 


Let V he as in Proposition 5.5| and fix t G [0,T). Then S{t) 

V 


Lemma 3.6. 

by (3.17), is self-adjoint and one can set Aq = (^(m,/ iq^, 
such that if X> Xq and v > vq, g & (0, qq), then one has 


Inl¬ 


and vq, — C'(nr, Aq 


, defined 


V 


(3.18) < (5(t)W,W)^, and |5(t)W|^, < Co|W|^,, 

where Cq^ = A^(1 + d), Cq = (^(m, Aq |V|^i „^). Moreover, S'{t) : ^ is well-defined and 

one has 


(3.19) vw G |5'(t)w|^, < c'(|(id - nf)aiV|^„„ + 0 |atV|^„„)|w|^,, 

where Cq = C'(m, Aq |V|^,,„). 

Proof. That S{t) is self-adjoint is obvious (recall in particular that Q(^) is orthogonal). Let us now 
decompose 

Q{D)S{Q{D))-^ = T"[Z] + i (Tst,,, + J + AA’^ 

+ ^ + Ts%^) - xiDy^^{T^[Z] - T“[0])x(n)i/2 

+ x(iA)i/2(T-[Z] - T"[0])x(n)i/2 _ (T-[Z] - T"[0]) 

= [Z] +5(1) + AAg ^ + 5(2) + 5(r), 

where Z is obtained from V by setting to zero all entries with n > A^ +1, and (notice Q(C)Z = Z) 


< 5 T(i)(t,x ,0 ='x(a(T“[Q(e)V(t,x)] -T“[Q( 0 Z]), 


<5T(2)(t,x,0 = x(a(T"[Z]-T“[0]). 


Using (3.14bl in Corollary 3.4 we can define Cq = (^(m, Aq , |Z|^„„) such that 

\f\N £ co\\N\l, < (T“[Z]W,W)^, and |T“[Z]W|^, < ico|W 


By (3.14e) in Corollary 3.4 and since V satisfies (3.11), one has for any (t, x, ^) G [0, T) xK'^xM'^\{0}, 
||T"^[Q(C)V(t,x)] - T"^[Q(C)Z(t,x)]|| < C(m, Aq \ |V(<,x)|)(i/"i + e|V(t,x)|). 


Similarly, the Lipschitz estimate (3.14e) yields for any |a| < 2 and (t,x,^) G [0,r) x x IR“, 

sup ||a|'dT(i)(t,x,,f)|| < ^(m, Aq \ |V|^^,a)(z/"i + £'|V|^^)(1 + |C|)"'“', 
xeR'^ 


so that the contribution from 5(i) is estimated thanks to Proposition A.3 (item i) as follows: 
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Similarly, using (3.14f) in Corollary 3.4 yields, uniformly in t G [Oj?’), 


|T^[Z]-T^[0]||^,.^ + ^||(T-[Z]-T"[0])nr| 




<C(V\tn,|V|_)|V| 




and one obt ains iden tical ly the corresponding estimates for derivatives with respect to Thus 

yield for any W G 


Propositions 


A.3 


and 


A.4 


|5(2)W|^, + -|5(2)n?^W|^, < C(m, |vL^)|v| 

Q 




W 


H~ 


One easily checks that the last contribution satisfies the same estimate: 

l‘5(r)W|^, + < C(m, |vL^)|v|^,.^|w|^_,. 

tr 


Here, we used that for scalar functions u G VP 


1,00 


) and w G H 


-Kmd 


), one has 


(3.20) 


< C'(d)|u| 






H-i’ 


where C{d) is a universal constant (by duality, since for any Lp G H^, one has vip G and 

{w,v(p)H-^-m < Cklij-ihlwi.-klj?!)' 

Thanks to the above estimates, one may choose A > Aq = 0 (111, ^ sufficiently 

large so that one has 

VWGL2(K<i)^(i+0, ((5(2)+5(,))W,W)^,+A(A-i(Z4)W,W)^, >0. 

It is now clear (using that Q(^) is orthogonal for ^ G \ {0}) that one can restrict g, v, A as in the 
statement so that (|3.18p holds. 


As for the second part of the statement, by definition of the quantization T (see Definition A.2) 
and since dt commutes with constant operators and Fourier multipliers, one has 

1 


Q(D)5'(t)(Q(D))-i = - (ra,(,T) +ro:(,T)) ■ 


Moreover, using (3.14e) in Corollary 3.4 one has 


|5t(T“[V(t,x)])|| < (^(m, Aq \ |V(t,x)|) (|(ld - nt)atV(t,x)| + p|i9tV(t,x)|) . 


Again, one easily obtains the corresponding estimates for derivatives with respect to ^ and Propo¬ 
sition 


A.3 (item i) yields 

ll'^ (^)||l 2-).L2 


< C(m,h^\ |V|^„) (|(ld - nf)9,V|^„ + g\dty\^^) 


Estimate (3.19) is proved, and the proof of Lemma 3.6 is complete. 


□ 


Lemma 3.7. Let V, W, R be as in Proposition\3.5\ and A, g, v as in Lemma\3.6[ Then one has 


(3.21) SdtW + {Q{D))-^ (^iiT"[0]B"[0]|D|+7)z^ Q(1?)W = 7^[V,W,R], 

where 

Z(t,x,e) = ix(6(T"[Q(OV(Cx)]B“[Q(e)V(t,x)]-T"[0]B"[0])|e|, 

and iR.[V, W, R] G satisfies 

| 7 Z[V,W,R]|^,(t) <Co|R|^,(t) + Ci|W |^,(0 

wit/i Co = C(m, Aq |V|^i ^), Cl = (^(m, Aq |V|^^) X Moreover, one has 

VtG [0,r), |5i(7)EW,W)^,| <C(m,AQi,|V|^„)|V|^,,^|W|',. 













18 


Vincent Duchene 


January 20, 2016 


Proof. Notice first that, using the rotational invariance property ( |3.9| ), one has 
B“[0]9,,W + B*'[0]9yW = i{Q{D)y^B^[0]Q{D)\D\y\l. 


Thus applying the operator S to (3.15) yields (3.21) with 


def 


(3.22) 7^[V,W,R] = 5R + {Q{D))-^Ti^Q{D)^ 
'1 


- mm- 


{TsT + VJ) + X^-:\D) Q{D) -B"[V]a,W+-B^[V]aj,W 


1 . 


1 , 


Q 


Q 


- {Q{D))-^T-[0]Q{D) { i(B"[V]-B"[0])5,W+i(B*'[V]-B^[0])5j,W) . 


We first note that, using estimate (3.14a) in Corollary 3.4 yields immediately 

1 , 


e 


A-^{D)Q{D) (B"[0]a,W+ B^[0]9yW)U = V^'(^)B"[0]Q(i^)l^|W| „ < C(m, V\A)|W 


0 


lL 2 - 


We then deduce from estimate (3.14c) in Corollary 3.4 that 

||B^[V] - B“[0]||^,_^ + IjB^^y] - B2'[0] ' 


live- 


< 0 C{hg ,m, |V|^^)|V| ; 


and in turn, by (3.20), 


vte [0,T), |A^i(n)Q(D) ('-B"[V]9,W+-B^[V]aj,W 

\0 0 

Thus there only remains to estimate 


lL 2 


<C'(m,VMv|_)|V| 




W 


i 2 - 


7^(l)[V, W] = {Q{D))-^r,r^^^Q{D)\N - ^{,Q{D))-^ {Tsj + T/y) Q{D) f iB"[V]9,W + iB«[V]9^W ) , 


\0 


0 


^( 2 )[V,W] (Q(iA))-i7-r,,)Q(i?)W 


- (Q(n))-iT^[ 0 ]Q(n) Q(B^[V] - B“[ 0 ])a,,W+ i(B^[V] - B^[0])9 j,W 


where 


Z(i)(t,x ,0 =' -<5T(t,x,0B"[Q(0V] and = -T"[0]5B(t,x,e). 

Q Q 


These terms are estimated exactly as in the proof of Lemma 3.6 i.e. using the paradifferential 
calculus of Propositions |A.3| and |A.4| together with the estimates of Corollary |3.4[ thus we do 
not detail. Let us just indicate why these contributions are uniformly bounded with respect to g 
small. The case of 7 ?.( 2 )[V, W] is quickly settled by (3.14c) in Corollary 3.4 As for [V, W], we 
decompose as above 

Ib^M = -B^[0] + Lb^M - B^fOl). 

0 0 0 

The contribution from the second component is uniformly bounded with respect to g small, again 
thanks to (3.14c) in Corollary 3.4 The contribution from the first component may also be uniformly 
bounded by remarking that 

||,5T(t,x,0B"[0]|| < ||(<5T(t,x,e)nf)B"[0]|| + ||5T(t,x,0((ld - nf)B"[0]) || < g x C{m,h^\ |V|)|V|, 


where we used ( 3.14a| ), ( |3.14e[ ) and (3.14f) in Corollary |3.4[ 
Altogether, one estimates TZ in (3.22) as desired, namely 

| 7 ^[V,W,R]L, <Co|RU+Ci|W 


lL2’ 


with Co = C(m, |V|^,.^),Ci = C{m,ho\ |V|^^) x |V 
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There remains to estimate 3?(7iEW, W)^2- By (3.14d) in Corollary 


3.4 


one has 


V(t,Oe[0,T)x 




vcc- 


<x(C)ie|C(m,V, V _) V 




and 5I(t,x,^) is symmetric. Proposition |A.3| (items ii. and hi.) as well as Proposition A.4 yield 

|2 


Vte [0,T), |3fi(7)zW,W)^2| <C^(m,V'.|vL„)|V| 

The proof of Lemma |3.7| is complete. 




W 


L2' 


□ 


Completion of the proof. Assume W is sufficiently regular, say W S C^([0, T); 
so that all the calculations below are well-defined. We compute the inner-product of iden¬ 
tity (3.21) with W; it follows 


^^('5W,W)^2<^|([9*,5]W,W)^2 


3?(z-(Q(i:'))-^T^[0]B“[0]Q(i:»)|D|W,W) 

Q 

3fi(7'rW,W)^2 + (7^[R,V,W],W) 


/L2 


The second term on the right-hand-side is identically zero since T“[0]B^[0] is symmetric. The 
other terms are estimated thanks to Cauchy-Schwarz inequality and Lemmata |3.6| and |3.7| Alto¬ 
gether, this shows that provided we restrict g G (0, po) and v > vq as in Lemma 3.6 one has 


(3.23) 


^^(5W,W)^2 < Ci|W|'2 + |mV,W]|^2|W|^2, 


with Cl = C'(m,/io \ V ) X ( V 




-I- p 9(V -I- (Id — nf)9tV ). Estimate (3.16) follows 


from Gronwall’s Lemma and the coercivity of S, i.e. (3.18) in Lemma 3.6 together with the 
control of 7^ G L^{Q,T\L?) provided in Lemma 3.7 The fact that the same estimate holds for 
general W G C°([0, T); L^(IR‘^))'^(^+‘^) solution to (3.15) may be obtained by density and thanks to 
a standard regularization process; see [IHl Theorem 7.1.11]. Proposition 3.5 is proved. 


3.3 energy estimate 


The energy estimate derived in Proposition 3.5 quickly induces a similar H‘^ estimate, for any 
s > d/2 by using once more the specific structure of our system of equations, namely that 
singular terms appear only as linear components of the system (13.41). 


Proposition 3.8. Let s > d/2 + 1 and \J,\N & C°([0,T); he such that V satis¬ 
fies (3.1 ),(3.11) with ho,v > Q and 9tV G L°°((0, T) x and 

dt\N + iB"^[V]a2:W -f iB2^[V]ayW = R, 

with R € A ssume that g G (0, po) o.nd v > vq with go,i^o o,s in Proposi¬ 

tion 3.5 Then one has, for any t G [0,T), 


(3.24) 


W|^,(t)<C'o(0)e' 


Cl y t 


W 


t=0 \l2 


Cl V (t-t': 


Co(t')|R|ff,(i') dt' 


with Co{t) = C(m,/io \ Ci = C(m,/io \ || and 


V 


def 




xR'i)’ 


Proof. Denote W® = AW G C'°([0, T); L2(Kd))Af(i+d)^ with A'* = (Id - A)®/^. Then one has 

(3.25) dt\l\P+ + -By[y]dy\N^ =A^R+ [A", -f [A", iB?'[V]]a,,W. 

g g ^ g ^ ^ g ^ 
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We shall apply Proposition [3]^ to the above system, thanks to the standard tools on Sobolev spaces 
recalled in Section [A.2[ Notice first that since since s > d/2 + 1, Sobolev embedding yields 


V 




< C V 




so that we only need to estimate the commutator on the right-hand-side of (|3.25|) to apply Propo¬ 
sition 


3.5 Since A® commutes with B“[0], one has 


[A",-B^[V]]a,,W = i[A",B^[V] - B“[0]]a,,W. 
Q Q 


Now, thanks to the product and commutator estimates recalled in Section rA.2[ and following the 
proof of Lemma [3^ one easily checks that, for any t e [0,T), 


|[A*,-B"[V]]a,W|^, = -|[A^B"[V]-B"[0]]a,W|^, <C(m,V\|V|^J|V|^,|a,W| 
Q Q 

Obviously, since B^^^] = Q((0, l))-iB^[Q((0,1))V]Q((0,1)) (by @ ), one has 


Hs-i- 


[A^ -B^[V]]9,W|^, < C(m, |V|^J|V|^,|9,W|^._,. 


One could apply the estimate of Proposition 3.5 to V,W®,A®R satisfying (3.25), but one 


obtains a slightly stronger result by stepping back and using directly the differential inequality of 
the proof, namely (3.23[): 


with Cl = C(m,/ig \ |V|^„) x -I- -|- |(ld - nf)atV|^,^), and 

|7^[R,V,W]|^, <Co|R|^,+a|W|^,, 

where Cg = C(m, fe g \ and |V||. 

Estimate (|3.24[) follows from Gronwall’s Lemma and the coercivity of 5, i.e. (3.18) inLemma3.6 


Again, this estimate is proved for sufficiently regular W® = A®W, but may be extended to general 
W® G T); solution to (3.25) by density and thanks to a standard regularization 

process. This concludes the proof of Proposition |3.8| □ 


4 Well-posedness and stability estimates 

In this section, we collect the information gathered in the previous sections, which quickly yield 
Theorem o as well as the propagation of well-prepared initial data (see Proposition |4.2[ ). 

Let us first recall that Proposition 2.1 offers the existence and uniqueness of a (maximal) solution 


to ( |1.4[ ), for sufficiently regular initial data. Our results give additional information on the large time 
behaviour of these solutions, using in particular the energy estimates of Proposition |3.5| a nd |3.8| 
These estimates, however, are based on a different formulation of the equations, namely ( |1.9[ ), which 
we claimed to be equivalent. Let us precisely state below in which sense. 


Lemma 4.1. Let U {g-\iX 2 , ■ ■ ■ Xn .uf,... € CO([0, T); il®(M‘^))^(i+'') 


with s > d/2 + 1, satisfying (3.1). Then\/^= {g ^( 1 X 2 , ■ ■ ■ ,Cn-,V 2 , ■ ■ ■ ,V 2 , ■ ■ ■ , 

defined by ( |1.7[ ), satisfies V G C"^([0, T); and 


V 


H‘ 


<Cim, VMULJIU 


\H^ 


Conversely, iJM G (^^([O, T); iL®(K'^))^*^^+'^) satisfies (3.1), then the change of variables (1.8) defines 
U G C'O([0,T);iL®(M‘^))^(i+®*) and 


|U|^, <C(m,V\|V|^J|V| 




Moreover, if D as above is a strong solution to (2.1) (or, equivalently, ( |1.4 [ ) ), then \l is a strong 
solution to (|3.4|) (or, equivalently, (1.9)^; and conversely. 
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Proof. Notice that U G C°{[0,T )-implies immediately V G C°([0, T); 

since is an algebra for any a > d/2. The converse is also tr ue si nce the multiplication by 

is continuous from 




to iJ" 


see Section 


A.2 


It follows, by continuous 


Sobolev embedding, that all the terms in (1.4) and (1.9) as well as in the calculations below are 
well-defined in C'°([0,T) x 

The evolution equations for v„ in (1.9) are straightforwardly deduced from the ones for u„ 
in (1.4). The evolution equation for w follows from 

N N 

= y] hidtUi + UidtiQi - C*+i) = hAn, - u^V • (hiU,). 


i=l 


i=l 


Plugging the expression for dtUi in (1.4), and using V • (hiUi 0 u^) = hi(ui ■ V)ui + u^V • (hiUi), 


yields immediately the desired expression of the evolution equation for w in ( |1.9[ ). 

The corresponding result concerning the compact matricial formulation of the systems, and in 
particular (3.5), is obvious by the chain rule. 

This proves the first part of the statement. The second part is identical since all these calculations 
are reversible: the Jacobian associated to the change of variables is invertible; see (3.3). □ 

4.1 Large time well-posedness; proof of Theorem |1.1| 


By Proposition 2.1, one can set t'o = ) such that if (1.6) holds with lyg, then there exi sts a 

unique U G C^([0,Tj„ax)', strong solution to (2.1) and U = U'". By Lemma |4 1 

the change of variable ([^ defines Vg CO ([0,r,nax);il' 


""(i+c*) and V satishes 


and V 


t=o 


dtM -f -B^[V]a,,V -h -B?'[V]5„V = 0, 
Q Q 


= C(m, hg j U‘" ) U“ . Let us denote 


T*(M) = sup{t G [0,r„iax), Ir ^ and ([L5|)-(fL6|) holds with /io/2, 2z/o}. 


We restrict our discussion below to M > |V 1^^,, so that (by continuity) T*{M) > 0. 

Using the system satisfied by V, (3.14a) in Corollary |3.4| as well as Sobolev embeddings, one 
checks 


Ivlll llvl 


L“(0,T;ff=) IK'ol ((0,T)xR-^) ^ ^'(m,/ig ,M)M, 


for any T G [0,T*(M)); recall Hf is defined in ( |3.6[ ). 
In particular, one has for any T G [0,T*(M)), 


N 


sup g\dt{g ^Ci)| + y'l^tCnl + |9t<| + < |||V|||^ < C(m,/io\M)M, 

(t,x)G[0,T]xR-^ 


from which we deduce 


VnG {1,...,A^}, \hn{t,-)-hn{0,-)\< [ |Cn+i - Cn| < C(m, V\x T 

and, similarly, 

VnG{2,...,iV}, \v/^\ + \vl\ < |<L^J + |u)(L^J+C(m,hg-\M)Mxr. 

It follows that, for given M, one can set vq, g//^ = C(m, M ) such that if G (0, go) and ( |1.6[ ) 
holds with then V satisfies the assumptions of Proposition 3.8 for any t G [0, min{T*(M), T**(M)}) 
with 


= C(m, ho\M)M. 
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We thus deduce the energy estimate 

Vt G [0,min{T*(M),r»(M)}), |V|^, < (Co(0)|U“|^, + 

with Co = C(m,and Cm = C{m,hQ^ ,M). In particular this shows that one can 
choose M = 2Co(0)|U“|^, such that 


< max {T-i,, C(m, U“) | U“|^.} • 


Going back to the original variables through (1.8) and by Lemma |4.1 
can restrict vq = C{m,hQ^, |U“|^J and the time interval [0,T), with T~^ 
that U is uniformly bounded and (2.3)-(2.4) remain satisfied. By continuity, uniqueness of the 


, this shows that one 
bounded as above, so 


maximal solution and thanks to the blow-up criteria in Proposition|2.1[ we deduce that one can set 


T~ 


<r-i <C(m,/roMU“Lj|U 


\H^' 


This concludes the proof of Theorem 0 


4.2 Propagation of well-prepared initial data 

Proposition 4.2. Let s > dl2 + 1 and U“ G satisfying ( |1.5[ )-( [L^ as in Theo- 

rem \l.l\ and denote U ■ ■ ■ XN,uf, ... G C°([0,T); 

the solution to (1.4) and U = U™. //U‘" satisfies initially 

+ V(g-'cnL. + • (W)L= ^ ^0, 


then there exists Co = C(Mo,m ,/iq such that 


n—1 


N 


VtG[0,r), |U|^,-f t|V(£» ^Ci)|i2 +-0 ^ < CoMoexp(CoMot), 


n—1 


uniformly with respect to g G (0, go)- 

Proof. Let u s de note by U G C°([0, T); the solution to (1.4) and U = U™ defined 

by Theorem 0 and by V G C°([0, T); the associated solution to ( |3.4[ ), namely 

dt\/ 4- -B“[V]a„V -f tB^[V]5yV = 0, 

Q Q 

obtained through the change of variables 0 (see Lemma [4.1[ ). Finally we denote W = dtd■ 
Using the above equation and the product estimates recalled in Section |A.2l one has immediately 
W G C C°([0,T) x R‘^). 

Differentiating the above system of equations yields 

-f tB"[V]a,W-G = -i5dB“[V])a,V - iat(B^[V])5„V. 

g g g g 


By construction (see the proof of Theorem 1.1 above), one can restrict £<q , i/q = C(m,/ig ,|U“|^J 
so that V satisfies the assumptions of Proposition 3.5 namely (3.1)-(3.11) for t G [0,r); and one 
has 


(4.1) 




By estimate (3.14a) in Corollary |3.4[ one has 
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What is more, the additional smallness assumption on V • yields 

uniformly for g G (0, po)- 

Finally, by (|3.14c ) in Corollary |3.4[ one has 


|a*(B"[V])a,V|^, + |9,(B«[V])9,V|^, < p |VL„o)|V|^ 

Altogether, after applying Proposition |3.5| with W = i9tV, one deduces 


5tV 


L 2 - 


(4.2) |9tV|^,(t) < CoMoe^‘>^°^ \l^+Co |5*V|^,(t') dt', 

with Co = C(m, /ig Mo). Applying Gronwall’s Lemma to |9tV|^2 (t) exp(— CoMof) yields 

Vt G [0,T), \dt\/\^2it) < C'oMoexp(C'oMot), 

for Co = C(m, /ig Mg) sufficiently large. 


Using again the system of equations satisfies by V, namely (1.9|, estimates (4.1) and (4.2) yield 
1 


-|e ^VCi|r2 + -|v ■ w| < C'oMoexp(C'oMot) 

g\ in g\ in 


for Cg = C(m, /ig , Mg) and uniformly with g G (0, go)- Proposition 4.2 is proved. 


□ 


5 Convergence results 

In this section, we investigate the asymptotic behaviour of the previously obtained solutions in the 


limit p —)■ 0. We first show that the solutions of the free-surface system (1.4) converge weakly 


towards solutions of the rigid-lid system (1.3). Strong convergence results are then obtained, first 


by assuming that the initial data is well-prepared, and then by approaching the oscillatory “defect” 
through rapidly propagating acoustic waves. 

5.1 Weak convergence: the rigid-lid limit 

Our first (weak) convergence result is the following. 

Proposition 5.1. As p —>■ 0, let U“ —>■ U“ G satisfying ( |1.5[ )-( [L6| ). Denote, for g 

sufficiently small, iC 2 ,g, • ■ •, Cv.g, Ui^g,..., G T); H'^the solution 

to (1.4) with C/g = U“. Then, as p —>■ 0, converges weakly (in the sense of distributions 

and up to a subsequence) towards ( 1 ^ 2 ,..., Cv: Ui,..., u^v)^ G L°°(0, T; H‘^a 

solution of the rigid-lid system (|1.3[), with initial data 


iRL I 


.0 = {( 2 ^ < - <5-inw“,..., u“ - 6-^UW^) 


where d Y)n=i^n, Y)n=iXri + C" - Of+iluif with convention C)" = Cn+i = 0 


n vA-iy- is the orthogonal projection onto irrotational vector fields. 

Proof. Restricting g G (0, pg) if necessary, U‘" satisfies the hypotheses of Theorem |l.l[ Thus one can 

define (C 2 ,e, • ■ • > Cv.g) Ui^^,..., Ujv,^)^ G (^^([O, T); from the solution to (1.4) 

with initial data U = U‘"; and 




with T ^,M = C{m,hQ , ) U 
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Thus, by Banach-Alaoglu theorem on L°°(0,T;= L^{0,T-,H ®(IR‘^))', we can extract 
a weakly converging subsequence (in the sense of distributions), that we still denote 1)^: 


U. 


def 


U = (C 2 , • • ■, Cjv, ui,..., ujv) 


with U e 

Let us first notice that since satisfies (1.41 one has, uniformly for g S (0, go), 
Vn G {2,..., A^}, ||^tCn,e||^oo(Q ^ C'(fri) ^0 ,M)M. 


Thus, since the embedding of in is locally compact and by Aubin-Lions lemma 

and Cantor’s diagonal argument, one has (again up to the extraction of a subsequence) ()n,g Cn 
strongly in (^^([O, T); It follows by the logarithmic convexity of Sobolev norms, that for any 

VnG {2,...,iV}, Cn,g^Cn in ([0,T);iLfJ . 

Thereafter, we fix s' G (1 + d/2, s) so that ^ (^^(IR) and is an algebra. For the same reasons 

as above, we have also 

Vn G {2,..., Af}, ^ u„ - u„_i in C°{[0,T); 

and, using (Id — n)VCi,e = 0, 

VnG {l,...,iV}, (ld-n)u„,g^(ld-n)u„ inC°([0,T);iVii)'' . 

Let us now define 


N 


N 


N 


^ ^ “1“ Cn,£» Cn+l,^) ^ ^ ^ ^ ^n—1,^) “1“ 


n—1 


n—1 


n=2 


Notice now that all the (quadratic) nonlinear terms are strongly convergent, so that 


N 


= y] ((5„ + Cn - Cn+l)u„ 

n—1 


and 


(Id - n)w, ^ (Id - n)w in c^{[o,Ty,H/J<^ 

with convention Ci = Ca^+i — 0- Notice also that passing to the (weak) limit in the first equation 
of (1.4) yields 

V • w = 0 and thus nw„ ^ 0. 


We now define 

Notice the identity 

N 


N 


— 6 ^Ilwg with S = S, 


N 


(d„ + Cn.e - Cn+I.(?)u„,e = Wg - y] hn,g S ^IlWg = (Id - n)wg - g{g ^Ci.e)<^ 


n—1 

so that 


N 


=' E - Cn+I.,)iin., ^ (Id - n)w in C°([0, T); hQ . 


It follows, using the formula ( [T^ with 7„u„,(,-7„_iu„_i_g = 7„(u„,g-u„_7g)+p2r„u„_i_g 

and Wg, and since all the components converge strongly in C^{[0,T); H/^y), that u„_g converges 
strongly as well. By virtue of uniqueness of the weak limit, one has 


u 


■n,g 


u„ inC°([0,T);iVf„J . 
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We conclude by plugging the decomposition + 5 ^IlWg into the evolution equations for 


velocities in (1.4). Notice the identity, using that S ^IIw^ is irrotational, 


(u„,g • V)u„_g = (u„_g • V)u„_g + (u„_g • V)u„,g - (u„_g • V)u„_g - 2^(1^ ^nw^p). 

Thus the only (quadratic) term which does not involve at least one strongly convergent factor turns 
out to be an exact gradient and independent of n; and so is the unbounded (linear) component of 
the equation, namely “V(p“^(('i^g). This shows, passing to the weak limit all the other terms in the 
equation, that there exists Vp, independent of n, such that 

n 

Vn e {1,..., A^}, dtUn + + (u„-V)u„ = -Vp. 

i^2 

It is straightforward to pass to the limit in the conservation of mass equations, so that 

N 

Vn e {2,... ,iV}, dtCn + ^V-(/iiUj) = 0, 


and we have already seen that V • w = 0. 

Thus drf (^ 2 ,..., Cw, ui,..., unV e L°°(0, T; H^) n C'°([0, T); hQ is a solution to (fOt, 


and one checks immediately that = ((^2 ",..., C]vi~^IIw™,..., — <5 . 


This concludes the proof of Proposition 5.1 


□ 


5.2 Strong convergence for well-prepared initial data 

Proposition 5.2. Using the notations of Proposition 5.1. let U“ —>■ U“ S as p —)■ 0, 

with U'" well-prepared as in Proposition 4-^ Then (up to extracting a subsequence) —?> 


strongly in C°([0, T); for all s' < s, as g ^ 0. 

Moreover, i9tu„,g ^ 5tU„ G and W(p“^Ci.e) 


VpHL g 


where Vp^'^ is the pressure associated with solution to (1.3). 


Proof. We follow the proof of Proposition 5.1 However, we may use additionally that, thanks to 


Proposition 4.2 (and using the system of equations (1.4) to control time derivatives) 


N 


\dt{Q ^Ci,e)||i2 + 


n—1 


It follows (up to the extraction of a subsequence) u„ and Hw^ —0 strongly in C°([0, T)-,Lf^^), 

and therefore in C°([0, T); for s' < s; and 9tU„,e ^ dtUn G L°°{0,T; L'^{R‘^)) (in the sense of 
distributions). 

By Banach-Alaoglu theorem, ^V(p“^Ci.e) ^ weak limit in L°°(0, T; L^(IR‘^)) when p —>■ 0, 

and passing to the weak limit in the velocity evolution equations shows that this limit is Vp^^. □ 

5.3 Strong convergence for ill-prepared initial data; proof of Theorem |1.2] 


The proof of Theorem |1.2| is divided in three parts. We first construct a “slow mode” approximate 
solution, thanks to the appropriate rigid-lid solution. We then construct a “fast mode” approximate 
solution, satisfying an acoustic wave equation with appropriate initial data. Finally, we show 
that, thanks to dispersive estimates on the fast mode, the coupling effects between the two modes 
are small, so that the superposition of the two components solves approximately the free-surface 


system (1.4) with appropriate initial data. The energy estimate of Proposition 3.5 applied to the 


difference between the exact and the approximate solution, allows to conclude. 
Construction of the slow mode. Using Proposition |5.2| with well-prepared initial data 


U“ =' (0, a",..., C, u“ - ..., uio - rinw“) 
















26 


Vincent Duchene 


January 20, 2016 


and artificially setting 0, one obtains in the limit 


a solution to (1.3) with initial data 

L.0 = (cr,..., ci^, u” - (5-inw“,..., - r inw“)^, 

where we recall that 6 is the total depth, w™ J2n=ii^ri + Cn ~ Cn’+i)u“ with 

convention = Cjv+i ~ ^ n VA“^V- the orthogonal projection onto irrotational vector 
fields. Moreover, the corresponding pressure satisfies G L°°(0,T; L^(R‘^)), and one has 


(5.1) 


U 


RLI 




<C(m,VMU“Lj|U 


Iff^ 


(since U^, the solutions of ( |1.4[ ) with = U“ from which is constructed, satisfy the same 

estimate by Theorem O- 

Let us prove that one may (uniquely) choose G L^(0,T; The regularity of 

strong solution to (|1.3[) is sufficient to claim that € L^(0, T; is uniquely defined by 


N 


N 


Sn + E ^ E 


where = <5™ + Cn~ Cn+i with convention = Cn+i = 0- Notice now that 

/ly (uy • v)uy + uy v • (h^uy) = v • (h^uy ® u^) 

and 

N n N N N / N 

y : tj" E n vcp=i: nvcy i: =i; c""+1: 


n—1 i—2 


2 = 2 


2=2 


It follows, using that is an algebra, that there exists g L°°{0,T; 

such that 

Ap^^ = + d^dyp^'y + dl^py^y. 

Thus one may define a unique solution p^^ g L°°(0, T; i/®(IR‘^)) by Fourier analysis. From the 
product estimates in Sobolev spaces (see Section A.2) and estimate (5.1), one has 

1 ^* lL“(0,T;H»(K‘i)) — ^0 i|U Iz/JlU |^„. 


We denote = (pp^^^ ^rl^ _ _ ^ ^rl^ ^rl^ ^ ^ ^rl) g lo°(o, T; Ifpis chosen 

sufficiently small, then satisfies (1.5) and therefore the change of variables ( |1.7[ ) defines 
satisfying (see Lemma 4.1) 


(5.2) 


V 


■slow I 




<C(m,VMU“|^J|U 


\H‘' 


It is easy to check, using = 1 — since V • ( Yl,n=i = Oj f^at 

satisfies (1.9) up to a remainder term denoted g L°°(0, T; iL®“^(IR‘^)); and 


(5.3) 


DSloW I 




<pC(m,VMU“Lj|U 


\H^' 


^Let us remark incenditally that such a choice of pressure with bounded energy is allowed thanks to the Boussinesq 
approximation applied to the rigid-lid system. Without the Boussinesq approximation, generic initial conditions will 
generate horizontal pressure imbalances, which in turn yield an apparently paradoxical evolution in time of the total 
horizontal momentum; see [IS]. One can check that the total horizontal momentum is preserved for the free-surface 
system as well as for the rigid-lid system with Boussinesq approximation. 
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Construction of the fast mode. We constructed above an approximate solution of ( |1.4[ ) (in the 
sense of consistency), but which does not fulfil the required initial condition. We correct this defect 
through an explicit “fast mode”. Define 0,..., 0, 0,..., 0, where 

and 


(5.4) 


with initial condition Q 


dt{Q-Xn + -V-w"^ = 0, 


+-V(e-icr) = 0, 

Q 


1 lt=0 


= Cl" and 


= IIw™, so that 


(5.5) 


V 


fast I 


V 


slow 


t = 0 '' |l2 


< C(m,VMU-Lj|U 


177=’ 


where V“ is defined from U‘" after the change of variables ([Hf- 

The above is an acoustic wave equation and is well understood. The following results may be 
found in [3] for instance. There exists a unique solution (p“^Ci“, w®'^) G It 

satisfies IIw'^'^ = for any t G M, and 


(5.6) 


Vt G 


1 Aao I _j_ _ I c 

f '’1 177 = ^ A I 177 = 


1/2 


I Dl l 77 = ^ A I l 77 = 


^^(m, /7 o-MU“LJ|U" 


1/2 


177 '’' 


What is more, since d = 2, one has the Strichartz estimates 


— 1 /-ac I 




1 llLj’(R;L'!(R‘i)) 


+ w" 


Il/(R;L‘i(H=*)) 


<Cpi/P(|Ci"|^.+C|w-|^,) 


where p, q are admissible, namely 2 < p,q < oo and h + ^ ^ — cr and ^ . Set for 

instance p = q = 6 and tr = ^. By a scaling argument and differentiating once the system, one has 

(5-7) l|V'‘^^lLe(o.T;ivce(«.)) < C(m, |^J|U“|. 

It follows that we control quadratic nonlinearities as 


I II < 11 \/' 

r ® '' llLi( 0 ,T; 77 i(R‘i)) - 


fast I 


L6(0,T;W^i’6(]Rd)) 

<01/6 C(m,Vi,|U“Lj, 


X V 


fast I 


lL6/5(o,T;H^i-3(R6i)) 


where we used Holder inequality, then Sobolev embedding and (5.6)-(5.7), with the restriction on 
the time interval, 

T-i <C(m,/io-i,|U“|^J|U”|^,. 

We deduce that satisfies the equations (1.9) up to a remainder R such that 

'li( 0.T;L2(R=*)) ® (ff—^O). 


(5.8) 


R" 


Control of the coupling terms. Let us denote V^pp = . Let us first check that V^'Pp is 


an approximate solution to (1.9), in the sense of consistency. From the above, we have that 

av"pp + iB''[v’"ppiaa.v’"pp + -B^[wppia„v’"pp = r"'°"' + + r™"p' 

g Q 


where R®*°'^ and R^^®* have been defined and estimated previously; and 
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By estimate (3.14c I in Corollary 3.4 one has 


V' 


slow I 


J|V^ 


fast 


I V‘ 


slow I 


IJ/1- 


We deduce from ( |5.2[ ), ( |5.6| ) and ( |5.7[ ), proceeding as for ( |5.8[ ), 

iRCOUplI 


(5.9) 


lLi(0.T;L2(B‘i)) 


0 (p^O). 


U| 


Let us now denote U € T); the strong solution to (1.4) with initial data 

= U'" as defined by Theorem 1.1 and V the corresponding solution to (1.9) defined by the 


change of variables ([T^. By Theorem and Lemma [4.1[ one has 


(5.10) 


V 




<C'(m,V\ U“LJ U 


\H‘' 


Proceeding as above, we find that the difference between the exact and the approximate solution, 
W ll'" V'^PP - V, satisfies 

5tW + ^B“[V]a,,W + iBS'[V]9yW = + R“"p' + R, 


with 


R ^(B^y’^pp] - B“[v])9,,v‘"pp + ^(B^^y’^pp] - B?'y])ayV^pp. 


Again, by estimate (3.14c[) in Corollary |3.4[ one has 


(5.11) 


lL2 


<C(m,VMV‘^PP|^,.^,|V|,„)|V"PP-V| 


L2- 


We now apply Proposition |3.5| with W = V^pp — V, and deduce 


|vapp_v|^^(t) < |v^PP-V|^yo)e^'lll'^lll*+ [ e^'lll^ll|(‘”*'^C'o(t')|R'‘°"'+R^''"VR“'^P‘+R|^yt') dt', 

Jo 


with Co{t) = C{m,hQ \ ^(t)), Ci = C{m,hQ \ ||V 


IL“(0,TxK‘*) 


)• 


Using the above control on the initial data (5.5) and remainder terms (5.3),(5.8),(5.9),(5.11) and 
Gronwall’s Lemma, we obtain 


sup |wpp-v|^yt) 

te[o,T) 


0 (p^O). 


By the logarithmic convexity of Sobolev norms and estimates (|5.2[),(5.6),(5.10), it follows 


sup IV^PP-Vl ,(t)^0 

te[o,T) 


(p —> 0) 


for any s' < s. 


This concludes the proof of Theorem |1.2[ with the exception of the uniqueness of the strong 


solution of the rigid-lid system (1.3). However, given two solutions 


U 


“f aRL aRL r: 

— \Pj ) S2,j ) ■ • ■ ) SAf.jJ “1., 




,u^L)^ G L~(0,T;LC(My)^(i+‘^) {j = 1,2) 


with same initial data, we may construct as above (i.e. multiplying the pressure with a p pref¬ 
actor) Uf7,U|i7 G L°°(0,r;HyM‘^))^(i+‘^) n Wt’;;,“((0,T) x K'^), two families of approximate 
solutions of the free-surface system ( |3.4[ ), for arbitrarily small p. Applying Proposition |3.5| with the 
equation satisfied by the difference between the two solutions (after the change of variable (O)), 
using (5.2), (5.3) and taking the limit p —>■ 0, yields 
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A Notations, functional setting and technical tools 

A.l Notations 

We denote by ^(Ai, A 2 ,...) a non-negative constant depending on the parameters Ai, A 2 ,... and 
whose dependence on Xj is always assumed to be non-decreasing. It may also depend without 
acknowledgment on the horizontal dimension, d; number of layers, N; Sobolev index at stake, s. 

Given X a topological vector space, X' denotes its continuous dual, endowed with the strong 
topology. 

Given H a Hilbert space and T : —> i? a continuous linear operator, we denote T* its adjoint. 


For 1 < p < 00 , we denote the standard Lebesgue spaces associated with the norm 

1 

\f\Lp dx^ < 00 . 


The space L°° = consists of all essentially bounded, Lebesgue-measurable functions / 

endowed with the norm 

I/Loo =^esssupxgK<i|/WI < 00 . 

For fc e N, we denote by = {/ s.t. VO < |a| < k, 9“/ € (endowed 

with its canonical norm) where we use the standard multi-index notation for a-differentiation. 

We denote by the space of continuous functions on with continuous derivatives 

up to the order k, endowed with the same norm. 

The real inner product of any functions /i and /2 in the Hilbert space is denoted by 

(/i , 72)^2 [ /i(x)/ 2 (x) dx. 

For any real constant s € M, denotes the Sobolev space of all tempered distribu¬ 
tions, / e such that |/|//» = < 00 , where A is the Fourier multiplier 

A‘'=iyid-A)i/2 = (id+|i?|2)i/2. 


For any function u = u{t, x) defined on [0, T) x with T > 0, and any of the previously defined 
functional spaces, X, we denote L°°{0,T;X) the space of functions such that u{t, •) is controlled in 
X, uniformly for t € [0,r), and use double bar symbol for the associated norm: 

ILILoo( 0.T;X) = eSSSUPtg[o_j.)|M(t,-)|^ < 00. 

For fc e N, C^([0,T);Ar) denotes the space of X-valued continuous functions on [0,r) with 
continuous derivatives up to the order fc. Finally, we denote, in order to avoid confusions, 

= fFy“((0,T) X ml = {/, such that /, V/, dtf € L“((0,T) x ML}. 


We denote by (•, •) and ['I the Euclidean inner product and norm on vector space M^ or and 
II'll the corresponding induced norm on the space of N-hy-N matrices with real entries 

(the choice of the norms has little significance). 

If the entries belong to a Banach algebra X {e.g. X = IF^’°°(M‘^) or X = iJ®(M‘^) with s > d/2), 
then we denote (•, •)^, and ll'Hx corresponding inner product, vector and matrix norms. 
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A.2 Product and commutator estimates 

We quickly recall the standard product, Schauder and Kato-Ponce estimates in a > (i/2. 

Proofs or references concerning the following results may be found for instance in [25]. 

Consider scalar functions f,g G By Sobolev embedding, one has f,g & C'°(K‘^) and 

It follows that the product is well-defined and continuous. Moreover, one has fg G and 

^ I /* I I 1? I ■ 


Moreover, if F S C^(M) with k G N, k > a and satisfies F(0) = 0, then F{f) G and 

A consequence of the above is of particular importance in our setting. Let g G be such 

that 1 + g > ho > 0. Then for any / G G and 


/ I 




(it suffices to remark = f — and apply the Schauder estimate with F a smooth function 

such that F{X) = if \X\ < I — hg and F{X) = 0 if \X\ > 1 — /io/2). Finally, we have the 
celebrated Kato-Ponce estimate for commutators: 


I [A'",/] 9x51^2 < C'|/|^j<,|9x5|jj<,-i < C\f\jj,\g 




A.3 Paradifferential calculus 

In this section, we recall some results concerning Bony’s paradifferential calculus. We follow the 
definition and most of the notations of [29| , although the latter is restricted to scalar functions and 
operators. The generalizations to (finite dimensional) vector spaces brings no additional difficulty 
since each operation may be reduced to a linear combination of entrywise scalar operations. 

Definition A.l (Symbols). Let k G N and m € N. We denote P™ the space of locally bounded 
functions A(x,^) : x —>■ AIjv(i+d)(C) which are with respect to ^ and such that for any 
a G N^, X 9|'A(x, ^) belongs to ^^’“(K^^) and there exists a constant Ca such that 

yf G < c„(i + |^|)™-i“i, 

where we use the standard multi-index notation for a-differentiation. For A G P™, we denote 
Mr(A;n)‘'A' sup sup ||(1 + |e|)l“l-™9j“A(x, ^. 

|a|<n{GR'^ 

Given a symbol A G P™, one can associate a suitable paradifferential operator. 

Definition A.2 (Paradifferential operators). Let A G P™. Then we define the paradifferential 
operator Ta, with symbol A, by 

(A.l) yu G TaU{x) = {2Tr)-‘^ [ df, 

JTSL'i 

where is defined with FxTfj^{ri,f) = ^f{r],f)FxA{g,f) where FxA{g,^) is the (component-by- 
component) Fourier transform of A{x,f,) with respect to x, and ip is an admissible cut-off in the 
sense of ]2fA Definition 5.1.4], and whose expression does not need to be precised. 

We now give some results used in this work. 
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Proposition A.3 ([IH]) (5.25) and Theorems 6.1.4, 6.2.4). Let m gN and assume d < 2. 
i. If A G r™, then for any s G M., Ta ■ —>■ is bounded, and 

\\Ta\\hs^h.-,^<C{s)M^{A-2). 

a. If A G r™ and B G T™ , then TaTb — Tab ■ -G -i) is bounded, and 

IITaTb < C'(s)Mr(A;3)Mr(B;3). 

Hi. If A G r™, then (Ta)* — Tj : if® —>■ is bounded, and 

< C(s)Mr(A;3). 

We will also make use of the particular cases of Fourier multipliers and paraproducts: 
Proposition A. 4 ([IHIj Theorems 5.1.15,5.2.8). Let A G F™ with m,k gN 

i. If A depends only on A = A(^) G C'°°(K'^), then Ta = A{D) where A{D) is the Fourier 
multiplier associated with A, i.e. 

vu e l 2 (]r^)^(i+'^), e 7^(0 = A(pju{0 = 

ii. If A depends only on x; A = A(x) G 1F^’°°(K‘^) with k gN, then 

\\TaU - AU\\^, < Cik)\\A\\^,_^ 

V|a| < k, WTaO'^U - Ad^^UW^, < Cik)\\A\\^,^^ \U\^,. 


B Eigenstructure of our system 


In this section, we give some information on the eigenstructure of the operators at stake in the 
multilayer shallow water model ( |1.9| ), namely B^[V] defined in (3.4); the translation in terms of 
the initial formulation (1.1), or A“[U], is immediate through the similarity transform (3.5). More 


precisely, we show that the above matrices are semisimple provided that each layer’s depth is positive 
and the shear velocities are sufficiently small. This provides the complete eigenstructure of the full 
symbol +^^B^[V] thanks to the rotational invariance property (see Lemma 3.2). We use 

this result in order to construct a symmetrizer of the system with the desired properties described 
in Section [3.1| (Lemma |3.3| ), but we believe that such information is of independent interest. 

Indeed, despite numerous works on the subject, the available information on the domain of 
hyperbolicity and eigenstructure of the multilayer shallow-water model is very sparse outside of the 
one or two-layer situation. The one-layer case is very classical and there is no need to discuss the 
subject here. Pioneer work on the two-layer case, in the limit p <C 1 and dimension d = 1 setting, 
include [ 35 ] for the free-surface case and [ 17 ] for the rigid-lid situation. We let the reader refer to [ 5 ] 
and references therein for the case of p > po > 0- Additional, numerical information may be found 
in [T3]. Sufficient criteria for the hyperbolicity of the bi-fluidic shallow water model in the general 
situation of dimension d = 2 and free-surface situation are provided in [I5l|32], while the rigid-lid 
setting is treated in mM- Starting from N = 3 layers, explicit results become out of reach except 
for very specific situations; see miisiiiii. 

In the general case of N layers, the author provides in m a non-explicit sufficient condition 
for strict hyperbolicity in the case of dimension d = 1, and for g > go > 0, using that the 
eigenproblem, in absence of shear velocities, reduces to a finite-difference analogue of the Sturm- 
Liouville problem (a similar tridiagonal reduction already appeared much earlier in the literature; 
see |5| and references therein). In |31j . Monjarret provides a sufficient criterion for hyperbolicity, 
for d = 1 or d = 2, by exhibiting the symmetrizer recalled in Section [^ Very precise information 
concerning the eigenstructure are also provided in an asymptotic limit which does not fit in our 
situation when N > 3, since it requires a sharp scale separation for densities between each layer; 
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see (211 (4.1)]. In this section, we build upon these works, by showing that the strategy of [H] 
extends to dimension d = 2 and arbitrarily small density contrast g] under hypotheses on the flow 
which are fully consistent (although, again, non-explicit) with the criteria given in |31j . 

Roughly speaking, we show that provided that shear velocities are not too large, the nonlinear 
evolution of the flow maintains N modes of propagation. One of them is the barotropic mode, and 
is responsible for the singular time oscillations of the system in the limit —>■ 0. In our scaling, the 
wave speed of the N — 1 baroclinic modes are uniformly bounded from above and from below, and 
remain isolated if the shear velocities are sufficiently small. 

Let us fix /iQ > 0 and denote 


V/io = {(i? ^Cii C 2 j • ■ ■, G ^-^(1+0 such that (3.11 holds with ho}, 

^ho = {(l?~^Cii C 2 j ■ ■ ■, Cn, 0,.. ■, 0)^ e such that (3.11 holds with ho}- 

We first remark that in the case of dimension d = 2, there are N trivial eigenvalues of 4B“[V]. 


Lemma B.l. Let V € V/i^. There are N linearly independent eigenvectors of 4B“[V], with corre¬ 


sponding eigenvalue gff m“(V), as given by (1.8 1 , associated with rank-one eigenprojections 


(B.l) 


P°[V]1^'(J^[V])-in2Ar+„J^[V], 


where n.„ is the orthogonal projection onto the variable. 


Proof. The result is straightforward when using the block formulation of A®[U], (2.2), and the 
similarity transformation (j^, i.e. B"^[V] = (J^[V])-iA"^[F(V)]J'f’[V]. □ 

In the following Lemmata, we give sufficient conditions on V € Yho allowing to complete the 
basis of eigenvectors with distinct and real associated eigenvalues. 

Lemma B.2. Let Z S Zhg. Then for any g > 0, lB^[Z] has 2N distinct, real, non-zero eigenvalues, 
p±n{'Z-) = ±Ai„(Z), n G {I,..., N}, with 


-/ri(Z) < • • • < -AtAr(Z) < 0 < ^iv(Z) < • • • < pi{Z). 


Moreover, one can set g^ ,Co = (^(m, hg , |Z|) such that if g G (0, go), then 


(B.2a) 

(B.2b) 


Co^<g\tii(Z)\<Co and Cg-i < |m„(Z)| < Cq, 
Ihn(Z) - Mn-l(Z)| > (n>2). 


The associated eigenprojections satisfy (recall the definition of Hf in (3.6)^ 


(B.2c) 

(B.2d) 

(B.2e) 


||P±i[Z]|| + ||P±„[Z]|| <C'(m,hg-MZ|) (n>2), 

|p±i[z](id - n^ll + ||(id - nnp±i[z]|| < e C(m,ho 1, |Z|), 

||P±„[Z]n^|| + ||n^P±„[Z]|| <^.C(m,hoMz|) (n>2). 


All these objects are smooth with respect to g G (0, go) and Z G l^hg- 


Proof. The eigenvalue problem concerning lB^[Z] is related by (3.5) to the one for lA^[Z], which 
given the simple block-structure exhibited in (2.2) may be reduced, for Z G to the eigenvalue 
problem for the following tridiagonal matrix 


(HR)-i = D(e,)D(f)-iA^D(h-i)D(7)AD(eg) 

(recall notations and identities in the proof of Proposition |2.1| ). Equivalently, we consider the 
eigenvalue problem for 


je D(f)i/ 2 (|_|R)-iQ(^)-i /2 ^ D(eg)D(f)-i/ 2 A^D(h-i)D( 7 )AD(f)-i/ 2 Q,'g^)^ 
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Since is a real, symmetric tridiagonal matrix with non-zero (positive) off-diagonals entries [43j . 
there exists Ai < • • • < Xn and (xi,..., x^r) an orthonormal basis of such that 

Vn e iV}, T®x„=A„x„. 

Since 7 „, r„, 0 > 0, one may check that det„(T®) > 0 where det„(T®) is the determinant of the 

n-by-n upper-left submatrix {i.e. leading principal minor) of T®, from which we deduce Ai > 0 . 
Now, let us consider 

T° = D(eo)D(f)-i/2ATD(/i-i)D(7)AD(f)-i/2D(eo), 


i.e. the matrix obtained from T® by setting to zero the first row’s and column’s entries. Using the 
above analysis on the leading principal minor of order 1, one obtains immediately that there exists 
0 = A° < A 2 < • • • < A^ and x°, such that T°x° = A„x° (n = 1,..., A^). The above eigenvalues 
and eigenvectors depend continuously on the parameters which are, by assumption, bounded in a 
compact set; in particular there exists Co = C(m, |Z|), for a given N, such that 

Co-'<AO<---<A^<Co and |A° - A^i| > Cq-i (n = 2,..., iV). 


Again by continuity, one can set Pq ^ ^0 |^|) ^^at if S (0, ^> 0 ), then A„ satisfy the 

above (replacing Cq with 2Co). Furthermore, we have (augmenting Cq if necessary) 




-1 


< Ai = det(T^) X (nOLsAn) < g^Co. 


Going back to the original problem, one deduces from ( |3.5[ ) and \2.2\ that for any n G N, 
T, ±An^^^ is an eigenvalue of AB’“[Z]. This concludes the proof of the first part of the 


/i±. 

statement with (|B.2a|),(|B.2b|). 


Notice now that the above analysis is not restricted to > 0. Indeed, the formula for above 
defines a (complex) tridiagonal matrix for any g £ and is single-valued, as least for g sufficiently 
small. Since the off-diagonal entries do not vanish, T^’, and therefore AB’'[Z] has 2N distinct and 
non-zero eigenvalues for any g ^ 0. Moreover, we are in the situation of [5^ Theorem II.1.10], 
namely is normal for a sequence gi 0 (simply restricting to ft G M), and therefore Xn(g) and 
the associated eigenprojection are holomorphic around g = 0. 

The extra information concerning the corresponding eigenprojections in the limit of vanishing 
g are obtained thanks to standard perturbation theory [521 Chapter II. 1]. Indeed, from (1.9) (see 


also the proof of Lemma 3.2 1 , one can write for any £> G C, 


B"[Z] = L + 0<5B[Z] 


with 


L = 


/ Oat j N j Oat 

_aN_^ I 0^ I _0_|y_ 

\ Oat I Oat ~[ Oat 


/-O- 


N = 


I 


V 


Otv-i 


0 


!0/ 


a 


N 

i=i 


and (5B[Z] is smooth with respect to Z and holomorphic with respect to g, and satisfies 

||5B[Z]|| <C(m,V\|Z|). 


It is obvious that L has only two non-zero eigenvalues: 

Lx^ = ±^/ay^, 

with (Id — rif )x^ = 0. 

Using that gg.±i (resp. ±-v/a) is a simple eigenvalue of B’^p] (resp. L), we introduce the 
Dunford-Taylor integral for the eigenprojection 


P±i[Z] 


-1 

27ri 


^(L - r,\d)-\-g6B[Z]{l - f drj, 
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where r±i is a positively oriented closed curve enclosing £»/i±i as well as ±-^0 , but excluding 
the other eigenvalues of B®[Z] and L, namely 0 and gfi±n {n > 2). One can restrict |e| < qo, 
^ |Z|) such that |^?|||(5B[Z]|| max,,gr±||(L —??ld)“^|| < 1/2, and therefore the series in 

the Dunford-Taylor integral is uniformly convergent. In particular, 


1 /* c 

P±i[Z] = — / {L - r]\d)-^ dij + (-(5B[Z](L-77ld)~i)'' dr?. 

Ztti Jyyi 

The first term on the right hand side is exactly the eigenprojection onto x^, and the series in the 
second term is immediately estimated. We deduce 


|P±i[Z]|| < Cih^\g, Izl), ||(ld - nr)P±i[Z]|| + ||P±i[Z](ld - n?^)|| < P Cih^\g, Izl), 


thus the first part of (|B.2c ) and (B.2d) is proved. 


One cannot directly use the same technique for the other eigenvalues, as they correspond to an 
exceptional point of B^[Z] = L + gSB[Z], namely 0 is an eigenvalue of L with algebraic multiplicity 
A^(l + d) — 2, and splits for g 0 in 2N — 2 distinct eigenvalues, gp±n{Z),n > 2 (completed with 
the (d— 1)1V linearly independent eigenvectors given in Lemma B.l which remain in the kernel). As 
a consequence, the Dunford-Taylor integral around this group of eigenvalues only yields a control 
of the total projection, namely the sum of the corresponding eigenprojections; or, if one integrates 
around a single eigenvalue, one cannot in general ensure that the series is convergent, even for g 
small. 

However, we have seen that -\n(g) is holomorphic in g and converges towards A„(0) > 0 as 
p —>■ 0, so that gp±ni'Z-) = ±g\n{g)~^^‘^ is holomorphic near = 0. This means [22l Theorem II.1.8] 
that the exceptional point is not a branch point, and therefore P±„[Z] is single-valued, but still may 
have a pole at p = 0. Using now that 0 is a semisimple eigenvalue of the unperturbed operator, L, 
and that we have shown that gp±niX), the eigenvalues of B®[Z] splitting from 0 are simple (as they 
are one-dimensional), we may use the so-called reduction process, and deduce [52 Theorem II.2.3] 
that the associated spectral projections, P±„[Z], are actually holomorphic at p = 0. 

This shows (B.2cl, and by continuity that P±„[Z] is a projection onto the kernel of L. Thus 
(P±„[Z] )nf = rif (P±„[Z] |g^o ) = Ojv(i-i-d), and ( |B.2e| ) follows. Lemma B.2 is proved. □ 


We now deduce from Lemma [B.2[ thanks to standard perturbation theory, the corresponding 
information on the eigenvalue problem for any V € 

Lemma B.3. Let V G V/j^. Then one can set z/ = (^(m, , |V|) such that if g G (0, go) and 

V satisfies additionally 

WnG {2,...,N}, |<|-P |w((| < 

then the matrix iB^[V] is diagonalizable. In addition to the N “trivial” eigenvectors described in 
Lemma B.l 
such that 


iB^[V] has 2N eigenvectors corresponding to distinct and real eigenvalues, pL±n(y), 

P-iiy) < ■ ■ ■ < M-Ar(V) < /X7v(V) < • • • < ^i(V). 

Moreover, there exists Cq = (^(m,Ivl) such that 


(B.3a) Cq ^ < g\fj,±i{\/)\ < Co and |At±n(V)| < Co (n > 2), 

(B.3b) |^_Ar(V) - ^Ar(V)| > Cg"^ and |/r±n(V) -/x±(„_i)(V)| > Cp"^ (n > 2). 


The associated spectral projections are smooth with respect to V S '^ho', and, for any Vi, V 2 G "Vho, 


(B.3c) 

||P±i[Vi] 

-P±l[V2]| 

1 <Ci e|Vi-V2|, 


(B.3d) 

||(P±„[Vi]-P 

±n[V2])n^ 

1 < Cl p|Vi — V 2 I (n > 2), 


(B.3e) 

||P±n[Vl] 

-P±«[V2]| 

< Ci(|(id-nf)(Vi-V 2 )| + e|Vi-V 2 |) 

(n > 2), 

where Ci 

= C(m,hoMVi| 

,|V2|). 
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Proof. We shall use perturbation arguments, starting from the knowledge that, thanks to Lemma [B.2[ 
the non-zero eigenvalues of B^[Z], for any Z £ are simple. In what follows, we denote Z £ 
the vector obtained from setting to zero all entries of V with n > -|- 1. Identifying (1.9) 

with (3.4), we may write (improving the description provided in Lemma 3.2) 

(B.4) B"[V] = B"[Z] -f gK^\d + p5Bf[V] -f e<5B,[V], 

where one can choose = K^{\/) (for example = (X^^i Si)~^w^) such that i5Bs contains (at 
first order in g) only contributions from shear velocities V 2 ,. ■ ■ ,vn, while ^Bf contains a leading- 
order contribution in but only on “fast variable” entries. More precisely, one has 


(B.5) 

(B. 6 ) 

(B.7) 


|(5B,[V]|| + ||<5Bf[V]|| < C(m,/ioMv|)|V|, 

||5B,[V]|| < C{m,h^\ \y\){,y-^ + g\y\), 
||5Bf[V](ld-n^)|| < ^^C(m,Ao^|V|)|V|. 


Since all the non-trivial eigenvalues of B^[Z] (and therefore the ones of B[Z] -|- gK^\6) are simple, 
we may use the Dunford-Taylor integral 


(B. 8 ) 


P±n[V] = 


-1 

27ri 


/ VR(77)(-e((5B, + ^Bf)R(r?))'= dry, 

'r±. k=o 


where r±„ is a positively oriented closed curve enclosing the eigenvalue p/r±„(Z) -|- gK^, but ex¬ 
cluding the other eigenvalues of B[Z] + gK^\d, and 
(B.9) 

N r. r, oO 


R{r]) = {B^[Z] -t gK^\d - ryld)-^ = ^ - 


P„4Z1 


P-n'[Z] 


+ 


P'JZl 


7 ^^ gfin' (Z) -I- gK^ - rj (Z) -|- gK^ - rj gK^ - rj ’ 


with fi±n(Z), P±„[Z], P° [Z] have been defined in Lemmata B.I and 


B.2 


Let us first consider the case n = 1. From Lemma 


B.2 


there exists Cq = C (m, /ig ^, | Z |) such 
that if n = 1, one can choose r±i as the circle of center eyr±i(Z) and of radius . Using (B.5) 
and (B.9), one can restrict g £ (0, po) with = C(m,/ig^, |V|) such that the series in (B.8) is 
immediately convergent and 

||P±i[V] - P±i[Z]|| < - ^ p"(sup||R(ry)||)'+”(||<5B,|| + ||<5Bf||)” < p C{xa,h^\ |V|)|V|, 

TT r _j_ 

n—l 

since P±i[Z] = ^ R(ry) dry is the first term of the series. Since P±i[Z], P±i[V] are rank-one, 
one has 

eyr±i(V) = tr(B"[V]P±i[V]) = yr/r±i(Z) + tr (B"[V](P±i[V]-P±i[Z])) + tr ((B"[V]-B"[Z])P±i[Z]), 


and the upper and lower bound on 0 |/r±i(V)| in (B.3a) follows from the ones on p|yr±i(Z)| given 
in (B.2a), and the above estimate with (3.10|) in Lemma 3.2 


We now turn to the case 2 < n < N. In this case, we set r±„ as the circle of center y?/i±i(Z) 
and of radius gCQ^ (with Cq as in ( |B.2b[ )), thus we may only ensure 

sup||R(ry)|| < e"iC(m, dg^ |Z|). 

r±n 


As a consequence, we need the precised estimates (B.6)-(B.7| as well as the decomposition into 
partial fraction (B.9) to ensure that the series converge. Indeed, thanks to (B. 6 ), one may augment 
1 / and lower pg in order to ensure 


g sup ||R(?y)||||(5Bs[V]|| < (^(m, dg \ |V|)(£)|V|-P rr ) < -. 


»?Gr±„ 
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Now, we use that for any n' > 2, 

jBf[v]p±„,[z] = (5Bf[v](id- n^))p±„,[z] + jBf[v](np±„,[z]). 


Thus by (B.2el and (B.7), one has 


|,5Bf[V]P±„.[Z]|| <C(m,V ,|V|) ^?|V| 


{n' > 2). 


The contribution from P®, [Z] is straightforwardly estimated, and the contribution from P±„' with 
n' = 1 contains no difficulty, since, (^3/i„(Z) + qK^ — ri)~^ is uniformly bounded. Altogether, one 
has 

sup ||e(5Bf[V]R(ry)|| < C(m, /iq ^, | V|) p|V|. 

»76r±„ 


Restricting g G (0, Qq) if necessary, the series in (B.8) converges and 

||P±n[V]-P±„[Z]|| < C(m,VMv|)(^?|V|+:/-!). 
In the same way, and using (B.2e|), one easily sees that 


and therefore 


pR(7;)nr <C(m,V\ V ) p V , 


(P±n[V]-P±„[Z])n?^|| < V )|V|. 


In particular, provided g G (0,^o) u ^ are sufficiently small, contains exactly one 
eigenvalue of iB^[V], and (B.3a)-(B.3b) follow. 

Estimates (B.3c),(B.3d),(B.3e) are obtained identically as above, using the decomposition 

B"[V2] = B"[Vi] + giK^ - K^)\d + .5B[Vi, V 2 ] 


where 


5B[Vi, V 2 ] =' B"[Z2] - B"[Zi] + e(<5Bf[V2] - 5Bf[Vi]) + e(<5B,[V2] - (5B,[Vi]), 

and remarking that 


\K- - K-\ + -||B"[Z 2 ] - B"[Zi]|| + ||5 Bs[V 2] - 5B,[Vi]|| < Ci(|(ld - nf)(Vi - V 2 )| + p|Vi - V 2 I), 

e 

|UBf[V2] - <5Bf[Vi]|| + -||(5Bf[V2] - ^Bf[Vi])(ld - nn|| < CAMi - V 2 I , 

Q 


with Cl = C(m, /ig |Vi I, IV 2 I). This concludes the proof of Lemma 


B.3 


□ 


Remark B.4. The proof of Lemma\B. 3| is somewhat cumbersome and rely on delicate properties of 


B^[V], namely (B.4) with (B.6) and ( |B.7 |, because we wanted to be as precise as possible as for the 
hyperbolicity conditions (see Remark 2.2). Th e pro of is considerably shortened and appears more 
robust if one replaces the assumption of Lemma B.3 with the more stringent |V| < v~^, as one may 
then simply use B'^y] = B^[Z] + p(5B[V] with |(5B[V]| < (^(m, hg |V|)|V| in lieu of (B.4). 
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